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ABSTRACT 


The  second  generation  design  of  the  Naval  Postgraduate 
School’s  ocean  floor  geomagnetic  data  collection  system  is 
described  with  emphasis  on  modifications  for  improvement 
over  its  predecessor.  These  improvements  include  15  channel 
parallel  analog- 1 o-serial  pulse  code  modulation  conversion, 
fiber  optic  data  link,  and  radio  telemetry  to  shore 
recording  equipment.  The  system  has  flexibility  in  operating 
depth,  an  increase  in  data  acquisition  time,  synchronization 
with  land  site  data,  and  a  data  farmat  that  is  readily 
converted  to  digital  for  computer  assisted  analysis. 
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i.  r  ntroductid^ 

This  thesis  research  is  part  of  a  continuing  effort  by 
the  Naval  Postgraduate  School  to  obtain  a  long  terra  data 
base  for  the  improved  understanding  and  interpretation  of 
the  ELF  electromagnetic  environment  in  the  vicinity  of  the 
sea  floor.  The  overall  project  emphasizes  the  importance  of 
acquiring  measurements  of  geomagnetic  fluctuations  on  the 
sea  floor  over  a  period  of  several  years  at  various 
locations  and  depths  while  simultaneously  obtaining  data  at 
a  land  based  site.  Primary  objectives  include  the 

interpretation  of  signals  covering  four  decades  of  frequency 
from  0.01HZ  to  100HZ  through  the  use  cf  total  field 
magnetometers,  induction  coils  and  OLF/ELF  receivers.  This 
phase  o,f  data  collection  covers  the  0.01PZ  to  20HZ  range 
utilizing  induction  coils  as  sensors. 

The  magnetic  noise  in  the  sea  is  of  interest  both  from  a 
geophysical  viewpoint  as  well  as  for  naval  applications.  In 
addition  to  undersea  field  magneton  at ry,  some  areas  of 
current  geophysical  interest  are  the  measurement  of  marine 
geomagnetic  anomalies  near  centers  cf  sea  floor  spreading, 
characterization  of  magnetic  fields  induced  by  ocean  waves, 
geophysical  exploration  and  the  utilization  of  low  frequency 
waveguides  present  under  the  sea  floor.  Applications  of 
interest  to  the  Navy  are  in  the  areas  of  mine-war f ar e, 
submarine  detection,  sab  marine  communication  using  a 
superconducting  quantum  interference  device  (S3(JID)  as  the 
sensing  element  for  an  extremely  low  frequency  (ELF) 
receiving  antenna  system,  and  possibly  as  a  means  of  remote 
verification  of  underground  nuclear  detonations. 


To  achieve  an  underst  and!  ng  of  the  nature  cf  the 
fluctuating  undersea  geomagnetic  spactrao,  it  is  necessary 
to  develop  a  data  base  with  goals  .or  high  quality  ana  great 
quantity.  The  collection  system  utilized  by  the  predecessors 
[ 1  ]  of  this  research  team  acknowledged  a  shortfall  in  both 
these  goals.  Their  collection  system  consisted  of  two 
chaanels,  each  with  a  coil  antenna  sensor,  orea mplif ier, 
variable  gain  amplifier,  voltage  controlled  oscillator, 
reference  oscillator,  and  an  analog  cassette  tape  recorder 
configured  as  depicted  in  Figure  1.1.  The  components  were 
encapsulated  in  Benthos  glass  spheres  with  che  deployed 
configuration  as  shown  in  figure  1.2.  The  data  collected 
with  this  system  was  limited  by  virtue  of  the  equipment 
restrictions  and  deployment  techniqaes.  The  analog  tape 
recorder  posed  the  greatest  limitation.  The  data  acquisition 
time  was  limited  by  the  cassette  tape  length  to  a  maximum  of 
45  minutes  per  channel  per  system  deployment.  The  accuracy 
of  the  data  was  greatly  dependent  upon  the  stability  of  the 
oscillators.  Additionally  the  utility  for  comparison  of  the 
data  with  data  taken  at  a  nearby  laud  station  is 
questionable  since  the  synchronization  of  the  data  could 
only  be  crudely  approximated. 

The  particular  objectives  cf  this  thesis  research  was  to 
improve  the  ocean  floor  geomagnetic  collection  system  design 
in  an  effort  to  improve  data  quality  and  quantity  while 
attaining  time  synchronization  of  land-cceaa  data.  A 
secondary  objective  was  to  develop  procedures  for  computer 
assisted  data  analysis.  Lastly  the  associated  thesis  report 
is  to  act  as  a  systems  manual  in  an  effort  to  provide 
continuity  for  future  system  users. 

In  the  design  of  an  ocean  floor  magnetic  data  collection 
system,  it  is  necessary  to  consider  the  constraints  placed 
upon  the  system  as  a  result  of  its  function  and  the 
environment  to  which  it  is  subjected. 


It  is  the  function  of  the  system  to  sense  geomagnetic 
fluctuations  near  the  ocean  floor  in  the  frequency  range  of 
interest,  process  and  transmit  tnose  signals  to  tno  surface 
for  telemetry  to  shore  based  receiving  equipment  for  dual 
(ocean-land)  recording,  decoding,  and  analysis  while 
maintaining  the  information  content  of  the  signals.  One 
requirement  of  the  integrated  systen  is  to  ensjra  signal 
purity.  That  is,  the  system  should  measure  the  magnetic 
signals  while  minimizing,  to  the  extent  possi'o system 
induced  variations  in  the  signal  being  measure!.  is  well 
known  that  a  magnetic  field  can  be  set  up  by  at  electric 
field  (Ampere's  Law).  This  mutual  interference  between 
system  components  must  oe  avoided  where  possible  and 
minimized  otherwise.  To  achieve  maximum  signal  purity,  the 
sensors  must  be  isolated  magnetically  as  much  as  possible 
from  the  rest  of  the  system.  3ecause  of  the  lew  signal 
strength  of  the  collected  data,  some  amplification  must  take 
place  at  the  sensor  location  to  allow  for  transmission  by 
nonferrous  metallic  cabLe  to  a  remotely  located 
instrumentation  subsystem,  further  magnetic  isolation  may  be 
achieved  by  utilization  of  a  fiber  optic  data  link  from  the 
instrumentation  sphere  to  the  telemetry  bouv.  Additionally, 
near  surface  interference  is  a  matter  of  concern.  Low 
frequency  electromagnetic  radiation  can  penetrate  the  ocean 
surface  for  up  to  several  meters  with  the  potential  of 
channelling  down  the  link  and  providing  perturbations  to  the 
siqnal  being  measured  if  a  metallic  conductor  cable  were 
used  for  data  transmission.  This  "antenna  effect"  can  be 
eliminated  by  employment  of  a  fiber  optic  link. 

Due  to  space  and  power  limitations,  transmission  line 
signal  attenuation  may  be  an  important  consideration  in  deep 
sea  locations.  Low  loss  fiber  optics  can  reduce  or  eliminate 
this  concern. 
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The  advantage  of  fiber-optic  links  have  be  an 
demonstrated  world-wide  by  experimental  optical  wave  guide 
systems.  Among  the  important  features  are  low  transmission 
loss,  wide  transmission  bandwidth,  insensitivity  to 
electromagnetic  and  radio  frequency  interference,  and 
suitability  for  digital  communications,  and  pulse  modulation 
methods  (for  amplitude  modulation  techniques,  fiber  optic 
cable  losses  are  independent  of  transmission  frequency) 
C21. 

Optical  fibers  are  attractive  since  a  single  strand  is 
sufficient  to  carry  the  information  signal  and  cannot  be 
short-circuited  by  exposure  to  sea  water.  Also  recent 
trends,  through  advanced  technical  expertise  in  fiber 
manufacturing  have  narrowed  the  economic  gap  between  long 
haul  optical  fibers  and  their  foremost  competitor,  copper 
wire  and  cables  [  3]. 

Ar.y  practical  design  ultimately  rests  upon  the 
requiremehts  of  the  user,  ais  resources  and  capabilities  as 
well  as  the  technology  of  the  day.  The  major  goals  have 
been  stated  as  that  of  signal  purity  and  data  quantity.  It 
has  been  shown,  although  not  rigorously,  that  a  fiber  op*  ic 
data  link  is  a  viable  alternative  to  the  metallic 
transmission  line  in  the  achievement  of  signal  purity.  The 
improvement  of  data  acquisition  time  is  gained  through  RF 
transmission  of  data  to  shore  receiving  equipment.  Data 
acquisition  time  is  thereby  limited  only  by  power  supply 
capacity  and  the  schedule  of  the  deploying  ship. 

As  indicated  earlier,  fiber  optic  systems  are  suitable 
for  digital  and  pulse  modulation  methods  of  data 
transmission.  In  addition,  analog  transmission  may  be 
accomplished.  There  are  several  pulse  modulation  methods 


available, 
pulse  width 


ues 


(PAM),  and  pulse  code  modulation  (PCS).  PCM  is 
known,  and  today  the  most  important  pulse  modulation  system. 
Three  separate  operations  are  involved  in  providing  a  PCM 
signal.  The  first  is  to  interrogate  the  message  signal  at 
regularly  spaced  intervals  (sampling).  The  second  is  to 
approximate  the  measured  amplitude  value  to  the  nearest, 
permitted  voltage  reference  level  (quantizing) .  The  third 
operation  is  to  represent  the  approximated  (quantized) 
amplitude  values  as  a  series  of  coded  pulses. 

In  PCR,  several  pulses  per  sampLa  are  used  to  signify 
the  amplitude  value,  instead  of  one  pulse  per  sample  as  in 
the  cases  of  PAM,  PPM,  or  PWM.  Consequently,  PCM  systems 
require  an  increased  bandwidth  (within  the  capability  of 
fiber  optic  systems)  for  their  transmission.  However,  any 
small  deformations  in  the  height  or  width  of  the  pulses  are 
irrelevant  since  it  is  only  necessary  to  know  whether  the 
pulse  is  present  or  absent  in  order  to  retrieve  the  origins' 
message.  Moreover,  in  a  PCM  transmission,  noise  is 
nonaccumuiat ive  because  noisy  PCM  signals  can  easily  os 
cleaned  up,  when  it  becomes  necessary,  by  the  process  of 
regeneration.  Thus  the  quality  of  a  PCM  transmission  is 
dependent  on  the  sampling,  quantizing,  and  coding  processes, 
and  not  the  length  nor  the  noise  of  the  transmission  media. 
By  contrast,  the  PAM ,  PWM,  and  PPM  systems  are  continuously 
affected  by  noise  and  cannot  be  cleaned  up,  or  regenerated. 
The  noise  is  accumulative  and  the  longer  the  distance  of 
tran smission ,  the  greater  will  be  the  noise  [4], 

Additionally,  the  desired  resolution  of  tie  measured 
signals  may  be  achieved  in  a  PCM  system  by  selection  of  the 
sampled  word  length  and  the  sampling  rate. 

Both  the  digital  and  analog  transmission  methods  are 
inferior  to  PCM.  The  loss  of  a  bit  of  digital  data  obscures 
the  value  of  the  sampled  data.  Analog  data  may  be 


18 


irretrievably  altered  by  noise.  PCS,  then,  is  perhaps  the 
ideal  transmission  method  for  our  application. 

This  report  delineates  an  ocean  floor  geomagnetic  dara 
collection  system  that  employs  coil  antenna  sensors,  a  pulse 
code  modulation  system,  fiber  optic  link,  RF  data  link,  and 
recording  equipment.  It  also  details  procedures  for  decoding 
and  digitizing  the  data  as  well  those  for  computer  assisted 
data  processing  for  eventual  analysis. 
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II.  3 1  ST  EM  DESIGN 


A.  EQUIPMENT  CONFIGURATION 

The  general  arrangement  of  the  Generation  II  System 
components  in  their  deployed  mode  is  shown  in  Figure  2.1. 

1.  Data  Collection  Eguj.£meiit 

A  fuctional  block  diagram  of  the  data  aquisition 
system  illustrated  in  Figure  2.2  reveals  the  following  major 
components: 

1.  coil  antenna  sensors  ( 2| 

2.  preamplifiers  (2| 

3.  preconditioners  (2) 

4.  pulse  code  modulation  system  (1) 

5.  optical  transmitt  er  (1) 

6.  optical  receiver  (1) 

7.  radio  frequency  transmitter  (1) 

8.  radio  frequency  receiver  (1| 

9.  instrumentation  tape  recorder  (1) 

The  system  components  are  discussed  in  the  following 
subs  eciior.s: 

a.  Sensor  Subsystem 

The  sensor  subsystem  is  currently  designed  with 
two  sensors  each  consisting  of  a  ceil  antenna  sensor,  a 
preamplifier  and  a  battery  power  supply  housed  in  a  Benthos 
glass  sphere  with  a  0.404  mater  inner  diameter.  The  sensor 
and  the  instrumentation  suosystsms  are  interconnected  by  two 
30  meter  long,  coaxial  cables  terminated  at  each  end  with 
Branrner  connectors  and  Bentos  sphere  penetrators  through 
the  glass  walls  to  the  enclosed  equipment. 
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stand  that  permits  ease  of  deployment  and  provides 
structural  support  and  protection  for  the  sensor  spheres. 
The  sensor  stand  design  calls  for  a  light,  free  flooding 
stricture  with  removable  strap  on  weights  to  obtain  the 
required  negative  buoyancy,  but  still  allow  for  portability 
and  ease  cf  handling  in  system  deployment.  A  aon-metallic 
material  is  necessary  in  the  stand  construction  to  avoid  the 
effects  of  induced  currents  due  to  electro-chemical 
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reactions  at  the  interface  of  a  metal  surface  and  the 
electrically  conductive  seawater.  The  stand  is  constructed 
cf  light  vsight,  high  strength. 
Aery lonitrils-Eut acl iene- St yrene  (A3SI  piping  which  is  well 
suited  for  this  application.  The  weights  are  fashioned  from 
lengths  cf  ASS  piping  filled  with  lead,  sand  and  transformer 
oil.  The  pipes  are  capped  off  at  each  end  to  provide  a  water 
tight,  r.cn-meta  Hi  o,  and  compact  means  of  adding  the 
necessary  negative  beoyarev  tc  the  stand. 

(2)  Coil  Antenna  Sensor  Each  sensor  is  a 
continuously  wound  coil  antenna  manufactured  from  5450  turns 
of  1 8  gauge  copper  magnet  wire  by  Elia  Engineering  of  Palo 
Alto,  California.  The  coils  weigh  a? pr oximately  100  pounds 
each  with  dimensions  as  depicted  in  Figure  2.4.  The  average 
enclosed  area  of  each  coil  is  3.0324  square  meters  which  is 
determined  by  taking  the  average  of  the  areas  obtained  by 
computing  the  area  lor  both  the  inner  radius  and  the  outer 
radius.  The  dimensions  cf  the  sensor  are  constrained  by  the 
geometry  cf  the  class  spheres.  The  coil  resistance  is  120 
ohms  with  a  self-indictar.ce  of  approximately  9.31  henries. 

(3)  Pree.mpli  fi er  The  preamplifier  selected  for 
use  in  the  sensor  sphere  is  the  model  13-10A  low  noise  ELF 
amplifier  manufactured  by  Dr.  ALan  Phillips  of  SHI 
International.  The  final  stage  of  the  amplifier  contains  an 
active  low-pass  filter  with  a  cutoff  frequency  at  20  Hz. 
Such  filtering  is  necessary  to  reject  the  effects  of  60  Hz 
signals  ar.d  harmonics  due  to  the  presence  of  power  lines  for 
the  similar  system  used  for  land  data  acquisition.  Stray 
frequencies  above  20  Hz  do  not  penetrate  signif lean- iy  to 
depths  a t  which  the  ocean  system  is  operated  (50-50  meters). 
The  overall  preamplifier  gain  is  approximately  55  15. 

A  schematic  diagram  of  the  preamplifier  circuit  is  contained 
in  Figure  A.  1.  The  gain  and  filtering  characteristics  of 


Figure  2.4  Induction  Coil  Sensor  Dimensions 


the  system  are  illustrated  in  Appendix  c.  and  ace  relevant 
to  the  discussion  concerning  system  transfer  functions.  The 
preamplifier  is  powered  by  two  .75  Ah',  12  volt  mercury 
batteries  mounted  in  a  section  of  1  inch  PVC.  The  batteries 
are  arranged  in  such  a  fasaion  as  to  give  the  required  plus 
or  sinus  12  volts  and  ground.  The  preamplifier  is  housed  in 
an  aluminum  case  to  shield  in  fields  developed  by  the 
preamplifier  circuitry.  The  preamplifier  and  battery  pact 
are  attached  to  either  side  of  the  coil  antenna  sensor  with 
electrical  tape.  The  interconnections  are  mads  by  coaxial 
connection  for  both  the  coil  input  to  the  preamplifier  and 
the  preamplifier  output  to  the  Beatnos  sphere  penetrators 
while  the  power  supply  input  is  made  with  a  keyed 
three-prong  Jones  connector. 

Power  drain  on  various  components  of  the  system 
is  a  serious  consideration  when  trying  to  maximize 


It  is  also  desirable  to 


on-statior.  data  collection  time, 
have  the  system  sealed  and  made  ready  for  sea  prior  to 
transportation  to  the  ship.  Thus,  an  externally  actuated 
magnetic  activation  circuit  is  employed.  This  circuit  shown 
in  Figure  A.  2  employs  a  normally  open  reed  switch.  The  reed 
switch  affixed  to  the  top  inside  of  the  glass  sphere  is 

closed  when  a  magnet  is  placed  in  its  vicinity  from  outside 

the  sphere.  The  closed  circuit  biases  transistor  switches 
closed  cutting  off  power  to  the  preamplifier.  This 

arrangement  is  advantageous  as  it  allows  packaging  ir. 

advance  of  the  of  the  deployment  date.  It  should  be  noted 
however,  that  this  circuit  does  deplete  a  nominal  amount  of 
power;  so  final  packaging  is  usually  completed  one  day  in 
advance  of  deployment. 


b.  Instrumentation  Subsystem 

The  function  of  the  instrumentation  subsystem  is 
to  receive  the  two  analog  outputs  from  the  preamplifiers, 
condition  them  to  the  requirements  of  the  V  ilow  on 
circuitry,  multiplex  the  channelized  signals,  convert  them 
from  analog  to  digital  format  and  pulse  code  modulate  them 
for  input  to  the  fiber  optic  data  lint. 

(1)  Equipment  Mo ua ting  Assembly  The  equipment 
mounting  assembly  is  shown  ir.  Figure  2.5.  Constructed  of 
alumninum  it  supports  two  Benthos  glass  spheres,  housed  in 
Benthos  hardhat  enclosures.  One  sphere  houses  the  instrumen¬ 
tation  subsystem  while  the  other  has  a  battery  pack  to  power 
the  optical  transmitter.  In  addition,  the  optical  trans¬ 
mitter  assembly  is  attached  to  the  vertical  member  of  the 
mounting  assembly.  The  design  incorporates  quick  release 
clamps  for  rapid  engagement  or  disengagement  from  the  main 
support  line  (250  feet  of  1  inch  polypropylene  line)  . 
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(2)  Signal  Conditioner  Pis  signal  conditioners 
receive  the  analog  signals  from  the  roil  preamplifiers  in 
parallel  fashion,  amplifies  them  on  ths  order  of  30  dB  and 
limits  signals  with  peak  amplitudes  greater  than  7.5  volts. 
The  circuit  is  presented  in  Figure  A.  3. 

(3)  Pulse  Code  Modulation  (PCK)  System  The  pulse 
cods  modulation  system  chosen  for  use  in  the  ocean  floor 
geomagnetic  data  collection  system  is  one  designed  and 
manufactured  by  Dr.  Robert  Lowe,  Lowecom  Inc.,  also  of 
Scripp's  Institute  of  Oceanography,  La  Jolla,  California. 
Tne  d.-vice,  schematically  illustrated  in  Figure  A.  4  and  A. 5 
features  15  channel  analog  input  capability  and  offers  the 
option  of  selectable  rates  of  2nsamples  per  second,  where  n 


may  take  on  integer  values  of  3  to  7.  By  appropriately 
jumpering  the  analog  input  pins,  the  sampling  rare  may  be 
increased  five  fold.  The  system  is  currently  deployed  with 
the  coil-x  output  connected  to  PCM  input  channels 
1,4,7,10,13  and  coil  y  output  connected  to  PCM  input 
channels  2,5,3,11,14  with  the  remainder  of  the  channels 
grounded.  With  a  selected  sampling  rate  of  32  samples  per 
second,  this  configuration  results  in  150  samples  per  second 
of  each  coil’s  signal.  While  channel  input  capacity  and 
sampling  rate  capability  are  much  greater  than  current  needs 
dictate,  they  provide  growth  potential.  The  additional 
channels  could  be  utilized  for  a  three  coil  system,  input 
from  a  pressure  sensor  for  swell  wave  data,  temperature 
sensor,  or  real  time  clock  data.  The  PCM  system  design 
includes  the  option  for  analog  or  digital  data  output  from 
channel  14  with  a  minor  change  in  the  external  pin 
interconnections. 

The  PCM  system  incorporates  a  crystal  oscillator 
and  associated  CMOS  integrated  circuitry  to  develop  the 
clocking  pulses,  a  16  channel  CIOS  analog  multiplexer,  a  15 
channel,  12  bit  CMOS  a nalog-to- 1 i g ita 1  converter  and 
associated  circuitry  to  provide  the  pulse  coding.  r he 
crystal  clock  oscillator  ooerating  at  a  frequency  cf  24.575 
KHz  produces  a  square  wave  output  vita  a  less  bit  rate  of  * 
bit  in  one  million.  This  frequency  combined  with  the 
inherent  delays  of  the  associated  clock-logic  circuitry 
provides  the  128  samples  per  second  sampling  cate.  By 
external  pin  interconnections,  the  other  sampling  rates  may 
be  selected.  The  clock  pulses  gate  the  analog  multiplexer, 
digi tal-to-analog  converter  and  follow  on  associated 
circuitry  that  form  the  pulse  code  words.  The  basic  output 
is  a  bi-phased  pulse  coded  signal  that  is  adjusted  to  TTL 
standards  by  adjustments  to  the  offset  and  gain  trimmers,  ?2 


and  P3  respectively.  The  data  is  organized  in  frames.  Eac r. 
frame  is  headed  by  a  sync  code  ward  which  is  followed 
sequentially  by  the  pulse  coded  samples  from  PCI  channels  I 
through  15.  The  sync  code  word  is  a  pulse  coded  digital  word 
with  a  decimal  value  between  0  and  4096.  This  word  is 
preselected  and  hardwired  on  the  circuit  board  by  connecting 
a  logic  high  or  lew,  as  appropriate,  to  pins  a  through  i  and 
k  through  m.  The  sync  code  word  is  essential  to  the  decoding 
process.  Additionally,  it  provides  protection  against 
confusing  the  land  acquired  data  with  that  taken  from  the 
ocean  flcor  as  each  system  has  its  own  sync  code  word. 
Currently  sync  code  words  35  38  ar.d  23  29  are  associated  with 
the  land  and  sea  sites  respectively.  3ync  code  word  3155  ha a 
been  used  in  test  and  sea  floor  data  collection  before  the 
September  runs.  The  3155  board  also  had  plus  or  minus  5  volt 
range  on  the  A/D  converter  while  all  other  boards  where  plus 
or  minus  10  volts. 

Gould  Gelyte  rechargable  batteries  (plus  and 
minus  12  Volt  and  grounii  are  utilized  for  all  power 
requirements  within  the  instrumentation  sphere.  The  current 
consumption  of  the  precond itioner  and  the  PC?!  system  is 
approximately  28  ma.  at  the  32  sample  per  second  rate. 

The  externally  actuated  magnetic  activation 
circuit  described  in  section  II.A.I.a.  (3)  is  utilized  to 
activate  the  instrumentation  sphere  equipment  upon 
deployment. 

The  equipment  interconnections  are  all  coaxial 
except  that  the  power  input  is  made  with  a  keyed  three 
pronged  plug  and  socket  connector.  Care  must  be  taken  to 
ensure  the  x  and  y  coil  output  leads  are  connected  to  their 
respective  PCM  inputs. 

Appendix  D  discusses  the  details  of  data  flow 
in  the  PCM  system  utilized  in  the  generation  II  geomagnetic 
data  collection  system  design. 
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c.  Optical  Data  Link 

The  optical  data  link  eliminates  the  concern  of 
antenna  effect  discussed  briefly  in  section  I  while  intro¬ 
ducing  a  variable  length  section  in  the  system  that  gives 
flexibility  for  operations  at  various  depths.  Fiber  optic 
system  design  considerations  as  well  as  an  overview  or.  the 
theory  ar.d  state  of. fiber  optic  technology  are  discussed  in 
Appendix  £. 
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Figure  2.6  Simplified  Block  Diagcai  of  the  Optical  Data 

Link 


The  fiber  optic  data  link  consists  of  an  optical 
transmitter  and  receiver  each  house!  in  an  environmental 
enclosure  and  connected  by  a  fiber  optic  cable  (see  Figure 
2.6)  .  During  deployment  the  cable  is  stored  on  and  dispensed 
from  a  plastic  reel  fashioned  froa  a  20  inch  bicycle  rim 
(see  figure  2.7  for  details!  . 
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(1)  Optical  Trinsmitrer/R?  reiver  Set  The  op-.:. 'a  1 
transmitter  and  receiver  are  the  3urr-Brown  3713T  and  ^ 7 1 .3 S. 
respectively.  The  3713T  and  3713R  when  connected  by  i 
suitable  fiber  optic  cable  form  a  250  Kbaud  NRZ  fiber  optic 
data  link  capable  of  operation  to  1.7  kilometer.  The  3713T 
fiber  optic  transmitter  is  an  elect: ical-t o-cpticai  trans¬ 
ducer  designed  for  digital  transmission  over  single  fiber 
channels.  Transmitter  circuitry  converts  TTL  level  inputs  to 
optical  pulses  at  data  rates  from  DZ  to  2  Kbaud  NRZ. 

The  3713R  fiber  optic  receiver  is  an  optical-* o- 
electrioai  transducer  designed  for  reception  of  digital  da-a 
over  single  fiber  channels.  The  receiver  circuitry  converts 
optical  pulses  to  TTL  level  outputs  with  a  receiver  sensi¬ 
tivity  of  15  nW  and  data  rates  to  250  Kbaud  NRZ. 

An  integrated  optical  connector  on  both  the 
3713T  and  3713R  allows  easy  interfacing  between  modules  and 
optical  fiber  without  problems  of  source/fiber/detector 
alignment.  The  metal  packages  of  the  3713T  and  3713R  provide 
immunity  to  electromagnetic  radiation  and  direct  prir.-sd 
circuit  board  mounting  with  no  additional  heat  sink 
required. 

A  simplified  block  diagram  of  the  3713T  trans¬ 
mitter  is  shown  in  Figure  A. 5.  Tie  input  stage  uses  a 
Schimitt  Trigger  Exclusive  OR  Gate  3,  for  noise  immunity  and 
its  logic  is  configured  so  the  phasing  of  .he  transmitter  is 
pin  programmable.  When  the  transmitter  phasing  terminal  is 
connected  to  theta-Q,  the  Light  output  is  in  phase  with  the 
digital  input  signal.  The  LED  is  on  when  the  TTL  input  is 
high.  Connecting  the  transmitter  phasing  terminal  to 
theta- 180  causes  the  reverse  to  happen:  the  LED  is  on  fcr  a 
digital  low.  This  option  is  selectable  on  the  ocean  floor 
geomagnetic  data  collection  system  tarough  the  installation 
of  a  single  pole,  double  throw  switch.  The  system  is 
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operated  in  the  in-phase  mode.  The  advantage  of  the 
theta- 180  mode  is  it  makes  it  possible  to  detect  a  Dreak  in 
the  fiber  cable  when  the  lata  link  is  idle.  However,  an 
idle  PCM  signal  does  not  provide  a  continuous  PTL  low  and 
thereby  obviates  the  advantage.  Amplifier  M  and  the 
current  switch  drive  a  light  emitting  diode  (LED).  The 
optical  power  output  of  the  3713T  transmitter  may  be 
adjusted  by  controlling  the  resistance  between  the  power 
adjust  pin  and  ground.  This  controls  the  peak  or  "on" 
current  in  the  LED.  When  the  resistor  is  minimum  the  LED 
output  is  maximum  and  as  the  resistance  is  increased,  the 
LED  output  asymptotically  approaches  a  minimum  value.  In  the 
ocean  system  the  resistor  is  adjusted  to  provide  a  240 
millivolt  peak-to-peak  receiver  output  signal. 

A  simplified  block  diagram  of  the  3713R  receiver 
is  shown  in  Figure  A. 7.  Input  light  is  converted  to  a 
current  by  the  PIN  photodiode  CR 1  which  is  connected  in  the 
photovoltaic  mode  for  maximum  sensitivity.  A  low  bias 
current  FET  input  currant- to-v cltage  converter  transforms 
the  diode  current  into  a  voltage  (7a)  which  is  further 
amplified  by  A2  and  presented  to  comparator  A3  as  Vb  where 
it  is  compared  to  the  threshold  voltage  Vt.  For  maximum 
noise  immunity  it  is  desirable  to  have  the  threshold  voltage 
set  to  a  value  corresponding  to  a  level  half  way  between  the 
higa  and  low  value  regardless  of  the  actual  light  level  at 
the  input.  In  the  3713R,  this  is  accomplished  by  a  peak 
detector  automatic  threshold  circuit.  A  pulse  of  light  input 
causes  a  voltage  pulse  at  7b  which  is  stored  ia  the  auto¬ 
matic  threshold  circuit,  divided  in  half,  and  supplied  to 
the  comparator  as  the  threshold  input  Vt.  Thus,  Vt  is  a 
voltage  corresponding  to  nne  midpoint  of  the  iighr  and  no 
light  conditions  of  the  diode.  Since  the  automatic  threshold 
cirouit  uses  a  capacitance  hold  technique,  the  threshold 


voltage  Vt  is  subject  to  decay  when  the  light  is  removed 
froa  CR 1 .  A  no-light  condition  of  approximately  0.5  second 
duration  (a  1  baud  data  rate)  can  be  used  with  no  significnt 
effect  on  noise  immunity.  The  PCM  system  presents  data  to 
the  optical  data  link  at  about  a  5K  baud  rata  and  is  there¬ 
fore  within  the  conditions  for  maximized  noise  immunity.  The 
analog  output  terminal  of  the  3713R  is  the  output  of  the 
linear  amplifier  A1.  The  voltage  at  ti i a  terminal  is  propor¬ 
tional  to  the  input  power  to  the  receiver.  As  such,  it  makes 
an  excellent  diagostic  point  for  testing  the  fiber  optic 
cable.  Monitoring  the  analog  output  terminal  gives  a  rela¬ 
tive  measure  of  cable  loss  at  the  transmitted  wavelength  and 
a  direct  measurement  of  receiver  sign al-to-noise  ratio  when 
the  transmitter  is  off.  In  the  ocean  floor  geomagnetic  data 
collection  system  the  modulation  input  is  taken  from  the 
optical  receiver  analog  output  since  the  output  amplitude 
can  be  adjusted  by  varying  the  optical  transmitter  power 
adjust  resistor.  An  output  level  of  233  to  240  millivolts  is 
utilized  to  ensure  RF  transmitter  modulation  without  over¬ 
driving  its  circuits. 

(2)  Fiber  Optic  Cable  The  fiber  optic  cable  is 
described  in  detail  in  appendix  E.  The  fiber  is  a  63.5 
micrometer  radius,  low  noise,  graded  ndex,  multimode, 
single  strand  optical  fiber  cladded  and  sheathed  with 
strength  members  as  illustrated  in  Figure  E.  2. 
Approximately  330  feet  of  type  5053  cable  is  stored  on  the 
storage  reel  shown  in  Figure  2. 7. a.  The  fiber  optic  link 
gives  flexibility  in  system  operating  depths  as  all  other 
system  parameters  are  fixed.  The  cable  is  affixed  to  the 
polypropylene  line  with  electrical  tape  at  about  ten  foot 
intervals  with  seme  slack  to  allow  for  st.reccni.ng  of  the 
main  support  line.  The  electrical  tape  is  fastened  in  such  a 
way  as  to  provide  a  tab  that  aids  in  quick  removal  during 
reco  very. 
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(3)  Couplings  The  coupling  system  ultiiized  is 
the  AMP  Optimate  Single  Position  Finer  Dptic  Cable  Connector 
System.  Since  AMP  connectors  were  not  available  with 


Figure  2.8  Piber  3ptic  Cable  Connector  Assembly 


parameters  compatible  with  the  fiber  optic  cable  intended 
for  che  system »s  optical  lata  link,  the  fiber  optic  cable 
was  giver,  an  additional  jacket  of  neat  shrink  tubing.  The 
heat  shrink  tubing  builds  up  the  outside  diameter  of  the 
cable  to  near  that  of  the  inside  diameter  of  the  connector 
as  shown  in  Figure2.8.  The  optic  cable  is  allowed  to  extend 
about  1/16  inch  from  the  face  of  the  connector.  The  cable 
and  the  heat  shrink  tubing  is  bonded  to  the  connector  shaft 
by  a  few  drops  of  cyanoacrylate  glue.  The  retaining  cap  is 
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secured  in  a  position  about  the  ferrule  by  the  crimp  ring 
mounted  1/16  inch  behind  tie  cap.  rha  polishing  bushing  is 
placed  ever  the  connector  and  securely  mated  to  the 
retaining  cap.  The  polishing  is  accomplished  by  applying  a 
figure  eight  motion  to  the  fiber  against  various  grit  abra¬ 
sive  paper  while  it  is  submerged  under  a  half  inch  of  water. 
A  9  X  9  inch  aluminum  baking  pan  is  very  suitable  for  the 
procedure.  Polishing  is  accomplished  using  5  grades  of 
silicon  carbide  finishing  paper;  3 23,  400,600,3  micron,  0.3 
micron  in  that  order.  The  polishing  is  completed  when  the 
surface  of  the  polishing  bushing  immediately  surrounding  the 
ferrule  is  polished  to  a  glossy  finish.  The  polishing 
bushing  should  not  be  reused. 

The  AMP  connectors  remove  concern  of  optical 
power  losses  due  to  lateral  misalignment,  angular  misalign¬ 
ment,  end  separation  and  end  preparation  quality. 

(4)  Bnvironmen  tal  Enclosure  The  optic  trans¬ 
mitter  and  receiver  are  each  housed  in  an  environmental 
enclosure  that  maintains  water-tight  integrity  in  the  ocean 
environment.  The  enclosure  illustrated  in  Figure  2.9  is 
manufactured  from  6  inch  I.D.  ?  VC  pipe.  One  end  of  the  pipe 
and  its  associated  end  cap  are  machined  to  a  true  round  as 
near  to  initial  dimensions  as  possible.  Typical  finished 
dimensions  are  ; 

PVC  PIPE  1.3.  6.09D  ±  0.003  inches 
0. D .  6.603  t  0.  003  inches 

PVC  END  CAP  I.D.  6.615  ±  0.015  inches 

An  "C'-ring  groove  is  machined  to  0.  39  3  inches  deep  by  0.153 
inches  wide  located  0.75  inches  from  the  machined  end  of  the 
pipe  to  conform  with  a  6.25  x  1/8  inch  "0"-ring.  The  pipe  is 
reinforced  with  rings  made  of  3.75  inch  sheet  PVC  material. 
These  rings  are  slotted  3.25  inch  below  ring  center  to  serve 
as  a  support  for  the  optical  transmitter  or  receiver  circuit 


boards.  The  machined  eni  cap  is  fitted  with  two  Brantner 
connectors  and  a  vent  plug.  The  opposite  end  cap  is  fitted 
with  the  fiber  optic  through-hull  penetrator.  Ail  these 
fittings  are  equipped  wita  an  "0"-riag  seal.  The  end  cap 
supporting  the  fiber  optic  through-hull  penetrator  is  bonded 
tc  the  unmachined  end  of  the  pipe  with  PVC  cement.  The  fiber 
optic  cable  is  passed  through  the  penotrator  and  initially 
fastened  with  cyanoacrylate  glue.  A  final  treatment  of 
DEV30N  5-minute  epoxy  serves  to  secure  the  cabla/penetrator 
interface. 

The  fiber  outside  the  penetrator  is  protected 
from  bending  stresses  by  application  of  three  layers  of  heat 
shrink  tubing  that  extend  Into  the  penetrator.  Each  layer  of 
tuning  is  longer  in  length  than  the  preceding  one  so  that 
the  cable  flexibility  increases  as  a  function  of  distance 
from  the  penetrator.  Initial  design  tests  indicated  the 
weakest  point  in  the  optical  link  was  the  point  at  which  the 
cable  entered  the  penetrator.  Tests  subsequent  to  the  appli¬ 
cation  of  the  tubing  indicated  this  point  to  be  one  of  the 
strongest. 

(5)  Battery  Poll!  Supplies  The  optic  transmitter 
is  powered  by  a  6  volt,  8.5  AH  rechargeable  sealed  lead  acid 
battery  configured  with  a  1  00  ohm  series  limiter  resistor. 
This  battery  is  boused  in  a  Benthos  glass  sphere  mounted 
with  the  instrumentation  sphere  in  the  equipment  mounting 
assembly  discussed  in  II .  A.  1  .b  .  ( 1)  . 

The  optic  receiver  power  supply  is  a  sould 
Gelyfe  (+  12,  -  12  volt  and  ground)  battery  pack  mounted  in 
the  can  buoy.  The  battery  pack  is  rated  for  18  A3 . 

d.  Can  Buoy 

The  can  buoy,  illustrated,  in  Figure  2.7, 
provides  71  pounds  of  positive  buoyancy  (this  is  inclusive 


of  battery  weight)  to  the  system.  It  contains  a  3ould  Gelyte 
rechargable  (♦  12,  -  12  volt  and  ground)  battery  pack  that 
powers  the  optic  receiver.  The  power  is  delivered  through  a 
three  conductor  cable  equipped  with  Brantner  connectors.  The 
car.  buoy  is  fitted  with  additional  3rantner  connectors  for 
data  transmission  in  and  oit  to  the  spar  buoy. 

The  axle  of  the  fiber  optic  storage  reel  is 
mounted  in  the  center  of  the  can  buoy  cover  containing  the 
electrical  connectors.  Tie  opposite  cover  is  fitted  with 
quick  release  clamps  which  are  used  to  secure  the  buoy  to 
the  main  support  line. 

e.  Telemetry  Spar  Buoy 


(1)  The  spit  buoy  designed  and 
constructed  by  Dr.  M .  3.  Thomas,  Applied  Physics  Laboratory, 
Johns  Hcpkir.s  University  and  LT .  P.  a.  Rutherford,  Naval 
Postgraduate  School  [5],  is  depicted  in  Figure  2.10.  The 
spar  buoy  is  partitioned  into  several  sections  as  follows: 

(a)  Transmitter  Power  Supply/Lead  Ballast 

Including  the  dampening  plate,  this  is 
the  bottom  two  feet  of  the  N  PS  spar  buoy.  It  is  constructed 
of  5  inch  PVC  tubing.  This  section  contains  three  6  volt/8.5 
AH  lead  acid  batteries,  which  provide  power  to  the  trans¬ 
mitter,  plus  ninety  pounds  of  lead  weights.  water  tight 
integrity  is  maintained  by  a  3/4  inch  thick  PVT  plug  with 
two  "Ow- rings  to  aid  in  pressure  sealing  (Figure  2.11).  The 
power  leads  from  the  batteries  are  soldered  to  a  pin 
assembly  permanently  housed  in  the  plug.  The  transmitter 
receives  power  through  a  cable  whin  attaches  to  the  pin 
assembly  in  the  plug  by  way  of  a  Brantner  connector.  The 
slip  coupling  adjoining  the  bottom  and  free  flooding 
sections  is  permanently  adhered  to  tie  bottom  section  with 
PVC  cement. 
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(b)  Free-Flooding  Section 

The  purpose  of  the  free  flooding  section 
is  primarily  to  enhance  the  separation  of  the  center  of 
buoyancy  (cb)  and  the  center  of  gravity  (eg).  This  section 
is  constructed  of  a  5  foot  length  of  5  inch  PVT  pipe.  k 
series  cf  holes  are  drilled  towards  the  bottom  of  the 
section  rc  act  as  flood  ports.  Holes  at  the  top  are  vents. 
The  power  cable  from  the  batteries  to  the  transmitter  tra¬ 
verses  this  section. 

( c)  Middle  Section 

This  section  is  composed  of  three 
sub-units,  the  lowest  portion  nousing  the  transmitter.  The 
transmitter  utilized  is  taken  froa  a  P3C  Orion  sonebuoy 
system.  It  has  a  power  output  of  one  watt  and  operates  in 
channel  29  at  a  frequency  of  171.65  MHz.  Its  watertight 
integrity  is  maintained  by  a  PVC  plug  identical  to  that 
between  the  lower  and  free  flooding  sections.  The  power 
cable  leading  from  the  bottom  plug  is  attached  to  this  plug 
with  a  Brantner  connector.  Precautionary  measures  that  are 
to  be  taken  during  this  portion  of  the  assembly  will  be 
discussed  in  the  next  section  entitled  transportation  and 
assembly.  Just  above  the  transmitter  housing  is  the  buoyancy 
collar  shown  in  Figure  2.12.  The  coLlar  provides  the  spar 
buoy  with  an  additional  eighty  pounds  of  positive  buoyancy, 
again  enhancing  further  separation  of  the  cb  and  eg.  The 
remaining  portion  of  the  niddle  section  is  composed  of  the 
middle  mast  which  houses  the  transmitting  cable  leading  to 
the  antenna.  This  middle  mast  is  reinforced  by  1/2  inch  and 
1  inch  PVC  tubing.  The  1/2  inch  sections  are  each 
individually  sealed.  The  1  inch  section  serves  as  a  channel 
for  the  transmitting  cable.  This  section  is  sealed  using  a 
PVC  ring  and  adhesive.  The  antenna  cable  protrudes  through 
the  center  and  the  hole  is  made  water  tight  with  application 
of  silicon  based  rubber  sealant. 
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Sealing  Plug  Handle 


(d)  Upper  Mast/A  rteana 

The  remaining  portion  of  the  telemetry 
spac  buoy  consists  of  an  sigh*  foot  section  of  1/2  inch  PVC 
tubing.  It  houses  the  remaining  length  of  the  transmission 
cable  and  is  capped  by  the  transmitting  antenna  which 
utilizes  a  reflecting  ground  plane  waich  is  at  150  degrees 
from  the  vertical. 

(2)  Trar.spqrta tic-g /Assembly  The  spar  buoy  is 
transferred  to  and  from  the  R/V  Acaaia  via  a  picitup  truck  in 
most  instances.  This  method  of  transportation  requires  that 
the  spar  buoy  be  easily  broken  down  into  lengths  that  are 
compatible  with  the  length  of  the  piccup  bed.  The  following 
subsections  discuss  the  procedures  for  the  assembly  of  each 
affected  section  and  the  precautions  required. 

(a)  Power/Sa  Hast  and  Free  Flooding  Sections 
Place  the  lead  ballast  in  the  bottom 

section  prior  to  inserting  the  batteries.  Ensure  that  the 
"Curings  are  sufficiently  coated  with  a  thin  film  of  high 
vacuum  grease.  Insert  plug  sc  that  it  snuggly  fits  below  the 
marked  level.  CAUTION:  prior  to  connecting  the  power  cable 
jack  to  the  plug,  ensure  that  the  indentation  in  the 
Brantner  plug  is  lined  up  with  its  counterpart  on  the  jack. 
Failure  to  do  so  could  damage  the  transmitter  by  reversing 
the  polarity.  Next  slip  tie  free  flooding  section  into  the 
joint.  The  holes  that  allow  flooding  are  to  be  at  the 
bottom.  Ensure  that  the  narks  or  both  slip  joint  and  free 
flood  section  are  aligned.  Bolt  the  free  flood  section  in 
place. 

( b)  Free  Flooding  and  Middle  Section 

Prior  to  joining  these  two  sections, 
insert  the  transmitter  assemoly.  First,  connect  the 
transmission  cable  then  slide  the  transmitter  in  place  again 
ensuring  a  snug  fit  and  proper  placement  of  the  plug  in  line 
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with  the  interior  markings.  CAUTION:  It  is  important  to  have 
a  53  ohm  lummy  load  or.  the  transmitter  at  this  point  of 
assembly  since  it  is  being  powerei  without  the  antenna 
attached. 

(c)  Middle  Section  to  Opper  Mast/Antenna 

This  portion  of  the  spar  buoy  is  to  be 
bolted  in  place  while  the  remaining  portion  is  hanging  over 
the  side  cf  the  ship  hoisted  up  by  the  Acar.ia's  on  board 
crane.  Ensure  that  the  transmission  cable  is  connected  prior 
to  upper  mast  attachment. 

(3)  Deployment /g ecovery /Pis assembly  All  sections 
of  the  spar  buoy  shall  be  assembled  onboard  the  Acacia  with 
exception  of  the  upper  mast  prior  to  deployment.  Attach  the 
crane  line  approximately  five  feet  up  from  the  tether  point 
on  the  spar  buoy.  Once  lowered  over  the  side,  the  buoy  can 
be  temporarily  secured  to  the  ship's  railing  and  lowered  an 
additional  amount  to  enable  attachment  of  the  upper  mast. 
The  buoy  is  next  lowers!  into  the  water  and  tended  along 
side  until  the  rest  of  the  system  components  are  deployed. 

Recovery  and  disassembly  is  accomplished  by 
following  the  reverse  of  tie  above  procedure. 

f.  Receiving  Station  Equipment 

The  receiving  station  equipment  consists  of  an 
AN/ARR-52A  radio  receiving  set,  an  H?  3964A  instrumentation 
tape  recorder  and  two  dual  trace  oscilloscopes.  AN/PRC  77 
radio  sets  serve  as  a  voice  communicat ions  link.  The  equip¬ 
ment  arrangement  is  as  shown  in  Figure  2.13.  The  antennas 
used  in  conjuction  with  the  telemerry  link  are  double,  tuned 
YAGI  type. 
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•  DATA  PROCESSING 

Currently  the  data  is  processed  to  develop  various  data 
plots  generated  on  the  Naval  Postgraduate  School's  IBM  3033 
computer.  Prior  to  the  actual  generation  of  the  plots 
however,  the  PCM  encoded  data  whim  was  stored  on  analog 
tape  utilizing  the  equipment  described  in  section  II.A.I.f. 
(receiving  station  equipment)  must  be  decoded,  digitalized, 
and  transcribed  onto  digital  tape. 

A.  HARDWARE 

1.  PCM  to  Digital  Convfrsion  Equipment 

Several  electronic  components  are  utilized  to  decode 
the  PCM  data  to  discrete  analog  form  and  then  to  digital 
fora  which  is  ultimatly  stored  on  9  track  digital  tape.  This 
data  processing  system  is  shown  schematically  in  Figure  3.  1, 
illustrating  data  flow  and  control  links.  Central  control 
of  this  process  is  accomplished  with  a  Hewlett  Packard 
9945A  computer  utilizing  an  operate  interactive  program 
titled  "PCM  PS3G".  After  execution  of  the  program,  the 
computer  requests  entry  of  specific  function  control  parame¬ 
ters  into  the  computer  and  other  equipment.  These  inputs  are 
used  to  control  synchronization  of  equipment  start,  digital 
tape  drive  speed,  decode  rate,  decode  time,  and  synchroniza¬ 
tion  code  word  entry  into  tne  decoder.  The  PCM  encoded  data 
is  fed  into  the  system  from  the  3P  3964A  4  channel  tape 
recorder  previously  used  to  store  the  data.  Utilization  of  a 
different  recorder  may  introduce  noise  due  to  a  difference 
in  head  alignment  and  drive  speeds.  Such  noise  decreases  the 
sian al-to-noise  ratio  upon  which  decode  synchronization 
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depends.  Should  decode  sync  hronization  not  occur,  a  lever 
analog  tape  drive  speed  and  sample  rate  must  be  selected, 
resulting  in  longer  processing  time.  The  decoding  of  the 
PCM  encoded  data  is  accomplished  with  a  Marine  Profiles, 
Inc.  Model  319  PCM  decoder.  Two  Monsanto  AM-64  1 9/OSM-368 
osci llosccpes  are  used  to  display  the  PCM  encoded  data  and 
the  discrete  analog  data.  The  discrete  analog  data  is 
processed  by  the  computer  to  develop  the  digital  data  which 
is  recorded  cn  6250  bpi/1230  ft.  9  track  digital  storage 
tape.  A  Kennedy  Model  9933  digital  recorder  and  computer 
interface  are  employed  for  that  purpose. 
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Figure  3.1  Decoding  System  Data  Flow  and  Control 


a.  Decoding  Procedure 


The  following  steps  illustrate  the  procedure  for 
decoding  the  PCM  encoded  raw  data  into  digitalized  data  and 
then  placing  it  onto  a  9  track  digital  tape  which  can  then 
be  delivered  to  the  IBM  3  03  3  computer  for  data  processing. 

1.  Disconnect  the  coaxial  cable  labled  number  7 
from  the  output  of  the  HP  39  68A  tape  recorder  located  in  the 
equipment  rack  adjacent  to  the  HP  9845A  computer. 

2.  Connect  cable  7  into  the  output  channel 

desired  for  decoding  of  the  HP  3954A  4  channel  tape 

recorder. 

3.  Energize  the  HP  3964A  tape  recorder,  the  two 
Monsanto  AM-64 19/USM-368  oscillo sopas,  the  Kennedy  interface 
and  AANDEEAA  tape  transfer  interface,  the  model  319  PCM 
decoder,  and  the  HP  9845A  computer  (the  ON/OFF  switch  is 
located  on  the  right  hand  side).  Also,  on  the  PCM  model  319 
decoder  place  the  fan  toggle  switch  in  the  OP  position  and 
the  AC/DC/OFF  toggle  switch  to  the  AC  position. 

4.  Place  into  the  right  hand  side  tape  reader  of 
the  HP  9645A  computer  the  program  named  "PCM  PR03". 

5.  Mount  the  analog  tape  onto  the  HP  3964A  tape 

recorder. 

6.  On  the  HP  9345A  computer,  type  the  command 
GET  "MT"  and  press  EXECUTE  to  initiate  processing. 

7.  On  the  PCM  decoder  place  the  following  func¬ 
tions  to  the  listed  positions: 

SOU  ECE  -  1 

SAMPLE  RATE  -  64  (for  3  3/4  recorder  speed) 

-  128  (for  7  1/2  recorder  speed) 

IN  V  ERT/NORM  AL  -  NORMAL 
OUTPUT/SAMPLE  RATE  -  3 
REC  CRDS/FIL  E  -  INFINITY 
SYNC  CODE  -  030 


8.  Press  RON  01  the  oomputer,  ignore  "enter  y  to 
skip  tape  ir.it"  {displayed  on  the  computer's  CRT)  ,  and  press 
CONTINUE. 

9.  The  computer  now  indicates  "load  tape"  into 
the  Kennedy  unit  and  "put  on  line".  To  do  this  energize  the 
Kennedy  ur.it,  load  the  rape  according  to  the  diagram  located 
on  the  inside  of  the  unit's  door,  press  the  LOAD  burton  and 
the  ON-LINE  button  located  on  the  front  on  the  Kennedy  unit. 

10.  The  computer  now  indicates  "enter  sync 
code".  Type  into  the  computer  3658  or  3155  as  appropriate 
for  the  sea  system  or  2620  for  the  Land  site  system  and 
press  CONTINUE. 

11.  Enter  transfer  time  in  minutes  and  seconds 
into  the  computer.  Example  20  minutes  and  20  seconds  would 
be  typed  in  as  20,20.  After  this  is  ions  press  CONTINUE. 

12.  Push  the  STOP  switch  on  the  PCM  unit  and 
push  CONTINUE  on  the  computer. 

13.  Push  the  START  switch  on  the  PCM  unit  simul¬ 
taneously  with  the  PLAY  button  on  the  Hewlett  Packard  tape 
recorder. 

14.  If  data  transfer  must  be  ended  early,  push 
the  KO  button  on  the  computer  to  haLr  data  transfer.  If 
this  option  is  elected  taen  T  nuet  be  entered  on  the 
computer  following  the  K0  command  to  write  end  of  tape  to 
the  rape. 

15.  "End  of  run"  will  be  indicated  on  the 
computer  CRT.  At  this  time  deenergize  all  the  equipment  and 
disconnect  cable  7  from  the  4  charnel  tape  recorder  and 
reconnect  it  to  the  8  chan  el  tape  recorder  located  in  the 
equipment  rack. 

16.  Finally,  remove  the  yellow  ring  from  the 
digital  tape  to  ensure  that  the  taps  is  nor  inadvertently 
over  written. 


2.  Main  Computer  Hardware 

"he  data  is  read  from  the  digital  tape  into  the  IBM 
3033  computer  via  software  developed  by  Lcdr.  J.  Fisher  of 
the  Naval  Postgraduate  School.  The  data  is  unpacked  (placed 
into  a  format  where  it  can  easily  be  accessed  and  manipu¬ 
lated)  ar.d  converted  into  a  spectrum  of  frequencies  using  a 
Fast  Fourier  Trc:nsform  routine.  This  data  is  then  plotted 
using  a  peripheral  Versater  plotter  for  spectral  display. 

B.  SOFTWARE 


The  programs  are  written  in  the  FORTRAN  IV  programming 
language  and  are  briefly  discussed  below. 

1  •  Test  Pro 2 ram 

The  test:  program  shown  in  Appendix  F  was  developed 
tc  test  the  computer  manipulation  against  known  data  to 
ensure  the  resulting  plots  would  meet  expectations.  With  the 
exception  of  the  internally  generated  data,  the  program  is 
essentially  the  seme  as  the  Main  Program  shown  in  Appendix  H 
and  explained  in  paragraph  3  below. 

2.  Mess  Storage  Program 

The  digital  data  tapes  are  mLabled  with  data  organ¬ 
ized  in  16  word  frames,  the  first  word  of  which  is  the  sync 
code  followed  by  15  data  words.  The  data  words  ace  organized 
x,y, z, x,y,z, .. ..  .  Since  the  sea  system  has  only  two  sensor 
coils  at  this  writting,  the  z  data  is  zero.  The  land  system 
data  has  actual  values  for  x,y,  and  z.  The  data  from  both  a 
land  and  sea  data  tape  is  eitered  into  the  IBM  3333  computer 
mass  storage  system  (MSS)  utilizing  the  mass  storage  program 
shown  in  Appendix  3.  Currently  only  tie  x  and  y  coil  data  is 
read  into  MSS  where  it  is  available  for  access  by  the  Main 
Program. 
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Main  Program 


The  Main  Program  utilizes  the  x  and  y  coil  inputs 
from  simultaneous  land  and  sea  data  to  develop  various  plots 
for  comparision  and  analysis.  The  Main  Program  utilizes  a 
large  number  of  arrays  with  the  intention  of  maintaining 
each  bit  of  original  data  or  calculation,  ready  for  recall 
at  any  time  for  father  computation.  Equivalent  arrays  are 
employed  because  the  plotting  routine  is  not  able  to  handle 

complex  arrays.  Ocean  arrays  are  indicated  with  the  suffix 

'O'  and  land  arrays  with  the  suffix  1  L'.  By  successive  calls 
of  11  SUBROUTINE  HD",  an  array  of  the  x  data  and  y  data  called 
xxx(I)  and  yyv  (I)  respect ively  is  established.  One  such 
array  is  called  a  block  of  data.  The  data  words  in  the  array 
are  integer  values  between  2048  and  4396. 

After  a  block  of  data  is  entered,  the  integer  values 
are  converted  to  complex  values  of  voltage  as  seen  at  the 
front  of  the  PCM  board  between  -13.3  and  +10.3  volts  for 
systems  utilizing  PCM  boards  with  sync  codes  of  2620  and 
3653.  For  systems  utilizing  sync  code  3155  the  voltages  vary 
between  -5.0  and  +5.0  volts.  The  xxx  and  yyy  arrays  repre¬ 
sent  data  taken  from  the  x  and  the  y  coils  respectively 

without  application  of  the  system  transfer  functions. 

The  voltage  values  (which  vary  over  time)  are  trans¬ 
formed  into  the  frequency  domain  by  utilizing  the  subroutine 
"FOUST"  which  performs  a  Past  Fourier  Transform  (FFT) .  It 
utilizes  the  Cooley-Tukey  FFT  algoritam. 

This  section  of  the  program  applies  the  system 
transfer  functions  to  the  data  that  has  been  transformed 
into  the  frequency  domain.  The  discussion  of  the  origin  of 
the  transfer  functions  is  given  in  Appendix  C.  The  actual 
transfer  functions  vary  as  a  function  of  frequency.  They  can 
however,  be  approximated  by  straight  line  segments 
different  frequency  ranges.  These  line  segments  are 
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generated  by  least  squares  approximation.  The  data  entering 
this  section  enters  with  units  of  voltage  and  is  converted 
into  units  of  nano-Teslas  oy  the  transfer  functions. 

The  previously  mentioned  sections  of  the  program  are 
embedded  within  a  do-loop  which  allows  the  analysis  of  a 
long  period  of  data  consisting  of  a  number  (NS)  of  blocks. 
The  Ith  element  of  each  block  is  taen  multiplied  by  its 
complex  conjugate,  averaged  and  converted  to  power 
(nano-Tesla  squared)  ar.a  stored.  Dnce  the  program  has 
passed  through  the  averaging  loop  for  the  last  time  it  has 
computed  the  arithmetic  average  for  each  frequency  point  in 
the  arrays. 

The  next  section  of  the  program  calculates  the 
Stokes  parameters  0  through  3  of  tie  individual  site  ortho- 
ganal  components  and  the  coherence  of  the  planar  circular 
polaiization  parameters  between  site  components.  The 
remainder  of  the  program  sets  ip  the  parameters  and  titles 
for  the  plots. 

A  more  detailed  description  of  these  programs  may  be 
found  in  Reference  [6]. 
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IV.  EXPERIMENTAL  RESULTS 

A.  DATA  QUALITY 

Test  data  points  were  developed  by  applying  a  5.0  volt 
peak-to-peak  sinusoidal  input  of  1,5,10,  and  15  Hz  to  the 
PCM  board.  The  data  was  sampled,  FCM  encoded,  transmitted, 
received,  taped,  decoded,  digitalized,  and  transcribed.  A 
computer  assisted  voltage  vs.  time  plot  of  the  1.0  Kz  data 
resulted  in  a  1.0  Hz  sinusoidal  signal.  A  plot  of  the  trans¬ 
form  of  the  5.0  volt,  5.0  Hz  signal  resulted  in  an  impulse 
function  at  5.0  Hz.  These  results  illustrate  than  the 
system  functions  as  designed  from  the  PCM  input  through  the 
data  processing  with  the  exception  of  the  application  of  the 
system  transfer  functions.  The  system  transfer  functions 
(which  take  into  account  tne  amplification  of  the  preampli¬ 
fier  and  siqnal  conditioners  as  well  as  the  transformation 
from  volts  to  nano-Teslas)  are  described  in  Appendix  C. 

3.  DATA  QUANTITY 

The  following  data  acquisition  period  resulted  in  the 
indicated  simultaneous  land/sea  data  from  La  Mesa  Village 
and  Monterey  Bay: 

1.  10:32  AM  to  2 : '4  3  PM,  25  Jily  1982 

2.  11:22  AM,  17  Aug.  1982  to  09:10  AM,  18  Aug.  1982 

3.  09:32  AM,  16  S=o.  1982  to  09:58  AM,  17  Sep.  1982 

This  represents  50  hours  and  30  minutes  of  accumulated  data 

during  three  collection  runs. 

Sections  A  and  B  above  illustrate  achievement  of  the 
stated  goal  to  develop  a  data  acquisition  system  with  a 
capability  to  gather  simultaneous  land/sea  data  with  high 
quality  in  great  quantity. 
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PLOTS  OF  PROCESSED  ACT'JA  L  DATA 


Appendix  I  consists  or  a  ouober  of  actual  data  plots.  It 
is  not  the  intent  here  to  analyse  toasa  plots  but,  to  merely 
illustrate  the  progress  of  the  Kav?.l  Postgraduate  School's 
Geomagnetic  Research  Team  effort.  On-going  work.  in  data 
analysis  includes  power  spectral  density,  coherence,  correl¬ 
ation  and  directionality  studies. 


V.  CONCLUSIONS  AND  RECOMM  ENDAMONS 


A.  CONCLUSIONS 

The  Generation  II  system  designed  for  this  research  was 
successfully  employed  in  the  continuing  effort  of  geomag¬ 
netic  data  collection  with  the  soacific  task  of  recording 
quality  simultaneous  land/sea  data  with  increased  quantity 
over  its  predecessor  systems.  The  data  collected  with  this 
system  was  only  limited  by  power  supply  capacity  and  by  the 
schedule  of  the  Research  Vessel  Acania. 

B.  RECOMMENDATIONS 

A  better  means  of  positioning  the  seniors  on  the 
seafloor  still  remains  a  problem.  The  ability  :o  accurately 
position  and  detect  sensor  orientation  will  -become  critical 
in  future  data  acquisition  and  component  field  correlation 
experiments.  The  exact  orientation  nethod  to  use  will  be 
directly  related  to  the  acquisition  system  chosen  for  future 
measurements.  A  flux  gate  magnetometer  may  possibly  be 
employed  to  indicate  orientation  of  the  antenna  coil  sensor 
assembly  to  the  earth's  magnetic  field. 

Althcuqh,  the  research  team  and  the  ship's  crew  were 
well  versed  in  the  system' s  configuration  and  deployment 
techniques,  the  deployment  time  averaged  about  43  minutes  as 
a  result  of  the  system's  size  and  complexity.  Furthermore, 
deployment  could  only  be  accomplished  in  calm  seas  at  a 
maximum  depth  of  200  feet.  To  alleviate  these  problems  the 
design  of  a  generation  III  system  is  recommended.  The  desian 
should  include  size  reducticn  and  system  simplification.  The 
simplification  may  be  achieved  by  concentration  of  several 
generation  III  components  into  a  single  package. 
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Additional  problems  to  overcome  in  future  systens  are 
the  lack  of  robustness  of  the  spar  o  u  o  y  and  filar  optic 
systems.  The  antenna  ground  plane  eLements  whece  manufac¬ 
tured  from  aluminum  and  ocassionally  fractured  under  stress 
during  transportation  or  deployment.  Future  systems  should 
utilize  stainless  steel  elements.  The  spar  buoy  sank  on  one 
occasion  due  to  compression  of  the  ouoyancy  collar  luring 
heavy  seas.  The  buoyancy  collar  should  be  replaced  by  a 
n one om pres sable  floatation  device.  During  recovery,  the 
fiber  optic  cable  was  crushed  on  one  occasion  -mi  kinked  on 
another.  This  problem  may  be  corrected  by  utilizing  a  fiber 
optic  cable  that  is  jacketed  with  a  heavy  rubber  s  heat. a. 

A  final  recommendation  is  to  incorporate  a  pressure 
sensor  to  measure  surface  and  internal  wave  activity.  This 
additional  information  will  allow  the  analyst,  to  determine 
the  contribution  by  waves  to  the  geomagnetic  data. 
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APPENDIX  A 


A.  EQUIPMENT  SCHEMATICS 

The  equipment  utilize!  in  this  project  * 
designed  and/or  constructed,  commercially  ava 
available  through  military  channels  and  locally 
The  system  equipment  schematics  are  present 
following  figures: 

1.  Preamplifier  (Figure  A.  It 

2.  Magnetic  Activation  Circuit  (Figure  A.  2) 

3.  Pre-amp  and  PCM  Preconditioner  (Figure  A. 3 

4.  pulse  Code  Modulation  System  (Figure  A. 4, 

5.  Optical  Transmitter  (Figure  A.  3) 

6.  Optical  Receiver  (Figure  A.7> 


Figure  A.  3  Pre-aip.  and  PCTfi  Preconditianer 
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APPENDIX  B 


A.  SYSTEM  AND  EQUIPMENT  PREPARATION 

In  order  for  equipment  to  perform  and  operate  as 
designed,  certain  preventive  maintenance  measures  and  proce¬ 
dural  steps  should  be  oaken.  Such  prudent  steps,  especially 
in  undersea  research,  will  most  certainly  allow  maximum 
results  and  prolong  equipment  life. 

The  following  list  serves  as  a  culmination  of  lessons- 
learned  in  previous  and  on  going  project  work: 

1.  Spherical  glass  instrument  spheres  are  composed  of 
two  Benthos  glass  hemispheres  which  have  precision 
ground  mating  surfaces.  As  ‘rhese  surfaces  are  susep- 
table  to  chipping,  care  must  he  taken  in  handling. 
During  assembly  of  tne  housing  prior  to  deployment, 
certain  procedures  are  necessary  to  maintain  watertight 
integrity: 

a.  The  glass  son  ere  should  be  dry  and  warmed  to 
room  temperature  using  a  heat  gur.  to  ensure  adhe¬ 
sion  of  the  sealing  tape. 

b.  All  mating  surfaces  must  be  cleaned  with  a 
solvent  such  as  alcohol  to  remove  any  foreign 
matter,  grease,  or  oil. 

c.  Prior  to  sealing  the  sphere,  careful  inspection 
of  all  equipment,  electrical  connections,  and 
packing  material  should  be  nade.  The  use  of  check 
off  sheets  and  coLor  coding  was  found  to  be  parti¬ 
cularly  helpful  in  avoiding  operator  error. 


2.  Electronic  equipment-  A  general  concern  is  the 
proper  application  of  power  supply  voltages.  Assembly 
of  a  large  number  of  electrical  components  within  the 
confines  of  a  small  glass  sphere,  necessitates  careful 
planning  and  utilization  of  United  power  resources. 
Improperly  applied  voltages  may  result  in  damage  to 
electronic  circuitry  and/or  loss  of  data.  Certain 
pieces  of  electronic  equipment  reruire  particular 
consideration  as  follows. 

a.  Sensor  system-  After  the  sensor  cools,  preamp¬ 
lifiers,  and  their  battery  power  supplies  are 
sealed  in  their  respective  spheres..  a  voltage 
measurement  should  be  t  a  ice  a  tc  ensure  operation 
and  then  the  deactivation  magnet  positioned  over 
the  reed  switch  such  that  the  voltage  drops  to 
zero. 

b.  _  PCM  system-  Extra  cars  should  be  taken  to 
ensure  the  x-coil  is  coni  acted  to  the  x-input 
channel  of  the  p  reco  r.aic i o  asr/a mpli f ier  ana  its 
x-cutput  is  connected  to  the  x-ir.pu:  channel  of 
the  PCM  system.  The  same  holds  for  the  y  and  z 
channels.  An  error  in  these  connections  will 
result  in  the  wrong  transfer  function  being 
applied  to  the  data  during  ta-o  computer  assisted 
data  analysis.  Tie  procedures  identified  for  the 
sensor  system  aoove,  hold  as  well  for  the  PCM 
system . 

c.  Tape  recorder-  The  recorder  heads  and  capstan 
drive  should  be  cleaned  aid  demagnetized  before 
each  acquisition  run.  This  action  will  enhance  the 
signal  to  noise  ratio  and  maximize  the  decoding 
speed.  High  quality  1300  foot  tapes  are 


recommended 


d.  Fiber  optic  lata  link-  Before  packaging  the 
fiber  optic  data  link  enclosure,  a  system  check 
should  be  run  to  determine  that  the  output  is 
20Q-240mv.  Before  connecting  the  power  cables  care 
should  be  taken  to  ensure  the  transmitter  power 
lead  has  the  ground  on  the  pin  marked  with  a 
dimple  on  the  connector  casing.  The  dimple  should 
be  aligned  with  the  indent  on  the  female  power 
connecter  mounted  in  the  end  cap  of  the  trans¬ 
mitter  enclosure.  Similar  steps  should  be  taken 
regarding  the  PCS  data  input.  This  concern  isn't 
necessary  with  the  receiver  power  input  as  it  is  a 
keyed  connector.  Once  the  connections  are  made 
remove  the  activating  magnet  from  the  PCM  sphere 
and  observe  the  idle  PCM  signal  on  an  oscillosope. 
Adjust  the  transmitter  power  adjust  potentiometer 
so  that  200-243  mv  is  observed.  Set  the 
0-degree/1 80-degre  e  phase  switch  to  :>.e  0-degree 
position.  Package  the  enclosures  as  follows: 

1)  Clean  "3 "-rings,  "0”-ring  grooves  and 
mating  surfaces  with  alcohol. 

2)  Lightly  lubricate  "0"-ring=  anc  pipe 
mating  surface  between  "0"-ring  and  pipe  end 
with  vacuum  grease. 

3)  Open  veut  plugs  an!  gently  slide  end  caps 

onto  the  enclosure  ensuring  the  electrical 
wires  are  clear.  CAUTION:  care  should  be 

taken  not  to  kink  the  fiber  optic  cable. 

4)  Both  enclosures  should  be  secured  to  the 
dispensing  reel  as  shown  in  Figure  B.2  and 
transported  as  a  separate  component. 
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3.  Batteries-  Gelyte  batteries  are  used  f sc  cost  power 
requirements.  These  batteries  should  be  checked  prior 
to  system  usage  and  a  systematic  approach  regarding 
service  use  and  charging  interval  should  be  taken.  The 
mercury  batteries  usei  in  the  sensor  spheres  provide  a 
constant  power  source  throughout  their  recommended 
life,  but  should  be  replaced  at  more  frequent  intervals 
to  ensure  successful  system  operation. 

4.  Connacto r/Cables  -  Before  and  after  each  use  in  a 


corrosive  environment. 


connectors 


should  be 


thoroughly  cleaned  to  extend  service  life.  Ail  rubber 
connectors  should  be  lightly  lubricated  with  silicon 
grease  to  ensure  electrical  continuity  and  prevent 
water  intrusion.  The  use  of  color  codir.a  on  cable  and 


Prevent 


connectors  lessen  the  possibility  of  error  and  is 
highly  recommended. 


B.  SYSTEM  DEPLOY  MENT/RECOVE  RY 

On  the  scheduled  day  of  system  deployment,  the  previ¬ 
ously  prepared  system  components  are  taken  to  the  Monterey 
Bay  Coast  Guard  Pier  by  pick  up  truck  where  trey  are  on 
loaded  to  the  R/V  Acania.  The  on  load  is  accomplished  with 
assistance  of  the  ship's  crew  and  its  on-board  crate. 
Figure  B.  1  shows  the  shipboard  layout  of  system  components 
that  has  proven  to  be  the  most  efficient  starting  point.  Not 
shown  in  the  figure  is  the  fiber  optic  data  link  and  storage 
reel.  It  is  temporarily  stored  in  tie  aft  deck  house  with 
its  interconnecting  electrical  cables. 

The  following  steps  are  taken  to  ? a cilit at e  final  prepa¬ 
ration  ar.d  deployment.  The  numbers  in  parenthesis  refer  to 
the  numbered  items  in  Figure  B.l. 
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Step  1:  The  crew  pisses  the  main  support  j.rae  *  1  > 
through  the  A-frame  boot  block  to  the  winch  on  top  or  the 
aft  deck  house.  The  line  is  reeled  on  tc  the  winch  drum 
allowing  enough  slack  so  the  line  may  easily  be  run  down  the 
port  side  of  the  ship  out  of  the  way  of  the  follow  on  opera¬ 
tions.  The  other  end  is  secured  to  the  sensor  system  (2) 
harness  by  means  of  a  bowline. 

Step  2:  The  spar  buoy  is  assamoied  starting  witr  the 
power  supply  section  (6)  and  the  free  flooding  section  . 
The  interconnecting  cable  (11)  is  ooanected  to  the  recep¬ 
tacle  in  the  power  supply  section's  water  tight  plug  and 
passed  through  the  free  flooding  section.  Remove  six  1/4-20 
screws  from  the  freef loodin g  section.  Align  the  mating  index 
marks  on  each  section  and  mate  the  sections.  Additional 
alignment  may  be  achieved  by  applying  a  long  saafted  screw 
driver  as  a  lever  in  the  free  flooding  ports.  Secure  t as 
sections  with  the  screws  previously  removed.  ,  Before 
connecting  the  power  cable  to  the  transmitter  section  (ft)., 
ensure  the  dummy  load  is  in  place  on  the  section's  antenna 
output  coaxial  cable.  Make  the  power  cable  connection.  Mate 
and  fasten  the  free  flooding  and  transmitter  sections  as 
before.  Do  not  assemble  the  antenna  mast  (9)  on  the.  spar 
buoy  at  this  point. 

Step  3:  Take  the  fiber  optic  transmitter,  receiver  and 
storage  reel,  packaged  as  shown  in  Figure  B.2  from  the  aft 
deck  house.  In  the  vicinity  of  the  cm  buoy  (5)  remove  the 
electrical  tape  securing  the  transmitter  to  the  storage 
reel.  Using  care  not  to  break  or  kink  the  optic  cable,  reel 
off  enough  cable  so  that  the  transmitter  may  be  positioned 
on  the  instrumentation  section  (4| .  Carefully  align  and 
thread  the  axle  into  the  hole  provided  in  the  canter  of  the 
electrical  connection  side  of  the  can  buoy.  Tighten  the 
fitting  until  it  is  snug.  Secure  tae  can  buoy  to  the  deck 
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(on  quick  release  clamp  side)  with  the  fiber  optic  cable 
protruding  from  the  power  supply  sphere  end  of  the  trans¬ 
mitter  enclosure.  Ensure  there  is  sufficient  room  for  elec¬ 
trical  cable  connections  at  the  PCM  sphere  and  cf  the 
enclosure  and  securely  fasten  the  transmitter  to  the  frame 
with  two  inch  Benthos  glass  sphere  sealing  tape. 


electrical  tape 


.r  ubbe  r 
wedge 


stainless  steel 
banding  strap 


Figure  B. 2 


Fiber  Optic  Enclosures  (1|  and  (2),  and  stora?' 
reel  (3)  . 


Step  4:  Connect  data  transmission  cable  (3)  to  toe 

sensor  subsection  (2)  using  color  code  for  easy  and  error 
free  connection.  Prior  to  making  connection  apply  a  thin 
coat  of  vacuum  grease  to  all  connector  shafts.  Make 
connections  and  secure  them  with  the  locking  sleeves.  Do  not 
screw  locking  sleeves  too  tightly  or  damage  to  the  sleeves 
may  occur.  Fasten  the  other  er.d  of  the  cable  to  the  cables 
protruding  from  the  pulse  code  modulation  (PCM)  sphere. 


Connect  the  PCM  data  cable  to  the  optic  transmitter  data 
input  connector.  Use  care  to  match  the  dimple  on  the  cable 
connecter  with  the  indent  on  the  connector  housed  in  the 
transmitter  enclosure.  This  is  essential  with  all  the  two 
conductor  Brantner  connections  in  tae  system.  Failure  to 
fellow  this  procedure  will  result  in  damage  to  the 
electrical  components  involved.  Connect  the  power  cable 
between  the  power  supply  sphere  and  the  optic  transmitter 
enclosure.  Now  connect  the  power  cable  between  the  car  buoy 
and  the  receiver  enclosure. 

Step  5:  Test  for  idle  PCM  signal  by  connecting  the  test 
cable  between  the  optic  receiver  enclosure  and  a:: 
oscilloscope.  Remove  the  power  actuating  magnet  on  the  top 
of  the  PCM  sphere.  Adjust  the  oscilloscope  to  observe  a  2 U C 
mv  PCM  signal.  The  signal  should  sync  reasonably  well. 
Remove  the  actuating  magnet  from  the  sensor  sphere  (one  at  a 
time)  and  observe  the  PCM  pulses  run.  The  system  is  now 
checked  out  from  sensors  through  the  fiber  optic  link. 
Remove  the  optic  receiver  data  output  test  cable  and  install 
the  data  cable  between  the  can  buoy  and  the  optic  receiver 
enclosure.  Install  the  spar  buoy  cable  (10)  to  the  can  buoy 
and  test  for  the  PCM  signal  by  placing  the  oscilloscope 
probes  on  pins  2  and  3.  Replace  the  magnets  on  the  sensor 
sphere  and  the  PCM  sphere.  Disconnect  the  spar  buoy  cable, 
the  optic  receiver  power  cable  and  tie  optic  receiver  data 
cable.  The  system  is  ready  for  deployment. 

Step  1  through  step  4  nay  be  accomplished  along  side  the 
pier  or  er.  route  to  the  deployment  iocation. 

Step  5:  Fasten  the  line  and  cable  to  the  spar  buoy.  When 
on  location  (180  foot  depta  for  this  example)  and  the  ship's 
anchor  is  properly  set,  the  ship's  crew  will  fasten  the  spar 
buoy  to  the  on-board  crane  and  hoist  it  over  the  starboard 
side.  When  the  fastener  is  at  the  rail,  a  timber  hitch  is 
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mads  about  the  spar  buoy  shaft  and  the  spar  buoy  is  made 
fast  to  the  ship's  rail.  Ths  crane  then  may  be  disconnected. 
Remove  the  four  1/4-20  hsithead  bolts  from  the  top  of  the 
spar  buoy.  Remove  the  dummy  load  from  the  end  of  the  coaxial 
cable  and  connect  the  cable  to  the  caole  in  the  antenna  mast 
section  (9).  Align  the  mating  index  marks  and  mate  the 
sections.  Secure  the  section  in  place  with  the  belts 
previously  removed.  The  rfrench  utilized  for  this  procedure 
should  be  fastened  to  the  operator  with  a  string  to  prevent 
loss  over  the  side.  When  the  sections  are  secure,  the  timber 
hitch  may  be  removed  and  the  spar  ouoy  lowered  into  the 
water  by  the  nylon  tending  line.  When  stable  in  the  water, 
remove  the  tape  on  the  spring  loaded  antenna  ground  plane. 
The  antenna  ground  plane  presents  an  aye  hazzard.  Use  cere 
in  tending  the  spar  buoy.  The  spa:  buoy  cable  should  be 
faired  along  the  starboard  side  outboard  of  any 
obstructions. 

Step  6:  secure  the  sensor  cable  connectors  to  one  of  the 
sensor  sphere  hardhats  using  electrical  tape.  Remove  tne 
actuating  magnets.  Allow  some  slack  and  tape  the  coaxial 
cable  (3)  to  the  main  support  line.  The  sensor  subsection  is 
hoisted  up  and  the  A-frame  adjusted  aft  until  the  subsection 
is  clear  of  the  ship's  stern.  The  sensor  subsection  is 
lowered  in  ten  foot  increments.  A  liberal  amount  of  slack 
must  be  allowed  between  fixed  points  along  the  coaxial 
cables  between  the  sensor  subsection  and  the  instrumentation 
subsection.  Electrical  tape  is  applied  at  ten  foot  intervals 
of  the  polypropylene  main  support  line.  A  tab  is  placed  on 
the  tape  to  allow  for  easy  removal  during  recovery. 

Step  7:  When  approximately  fifteen  feet  of  the  cable 
remain,  the  A-frame  boom  is  adjusted  in  so  that  the  main 
support  line  is  very  close  to  the  stern  of  the  ship.  An 
ample  amount  of  fiber  optic  cable  is  reeled  off  the  storage 


reel  to  allow  the  instrumentation  subsection  to  be 
positioned  on  the  stern  wita  PCS  end  down  and  quick  release 
clamps  toward  the  line.  The  main  support  line  is  lowered  to 
where  there  is  sufficient  slack  in  the  last  section  of  the 
coaxial  cable.  The  cable  connectors  are  fastened  to  the 
Instrumentation  subsection  frame  with  electrical  tape  and 
the  actuating  magnet  is  removed  from  the  PCM  sphere.  The 
quick  release  clamps  are  made  fast  to  the  main  support  line. 
The  main  support  line  is  hoisted  off  the  deck  and  again  the 
A-frame  boom  is  adjusted  aft  until  tie  load  is  clear  of  the 
stern.  The  load  is  lowered  again  in  ten  foot  increments  and 
the  optic  cable  is  fastened  to  the  main  support  line  with 
electrical  tape  as  before.  Care  must  be  taken  to  prevent 
breaks  or  kinks  in  the  fiber  optic  cable.  Coordination 
between  lowering  the  load  and  reeling  off  the  fiber  optic  is 
a  must.  When  the  170  foot  marker  (two  black  stripes  (150 
ft. |  followed  by  two  successive  single  black  stripes  or.  the 
main  support  line)  is  at  the  A-frame  boom  block,  the  winch 
lowering  the  system  is  stopped.  The  remaining  fiber  optic 
cable  is  secured  to  the  storage  reel  with  a  rubber  wedge 
inserted  in  the  storage  groove  and  taped  in  place  with 
electrical  tape.  The  storage  real  is  immobilized  by 
fastening  a  piece  of  parachute  cord  from  the  can  buoy 
strongback  through  the  storage  reel  spokes.  The  power  and 
data  cable  connections  are  made  between  the  can  ouoy  and  the 
optic  receiver  enclosure.  The  spar  buoy  tending  line  is 
fastened  to  the  can  buoy  by  means  of  a  shackle  and  the  spar 
buoy  cable  is  connected  to  the  can  buoy.  The  main  support 
line  is  lowered  until  the  170  foot  nark  is  at  deck  level. 
The  excess  fiber  optic  cable  must  be  carefully  tended  at 
this  point.  The  can  buoy  quick  release  clamps  ace  made  fast 
to  the  main  support  line.  The  Load  is  hoisted  and  the 
A-frame  boom  adjusted  aft  so  the  load  is  clear  of  the  stern. 
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The  can  buoy  is  lowered  into  the  water  with  the  excess  fiber 
optic  cable.  Simultaneously  the  spar  buoy  cable  is  released 
over  the  side.  The  end  of  the  polypropylene  main  support 
line  is  unfastened  from  the  winch  cable  and  fastened  to  a 
marker  buoy  then  released  over  the  stern.  The  system  is  now 
deployed. 

Step  8:  Check  with  tie  receiving  station  personnel  to 
ensure  they  are  receiving  a  sea  side  ?CM  signal. 

Step  9:  The  recovery  operation  is  just  the  reverse  of 
the  deployment  operation  omitting  the  system  checkout 
procedures. 
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APPENDIX  C 


A.  SYSTEM  TRANSFER  FUNCIIDNS 
1  •  Introduc t  ion 

The  triaxial  coil  magnetometer  system  is  diagramed 
in  Figure  C.  1.  Presently  only  the  Land  based  system  employs 
three  coils,  the  ocean  system  uses  two  coils.  The  data 
sensed  by  this  ma gnetcmet ec  system  is  given  by  Vx(t),  Vy(t) 
and  V z(t),  the  output  voltages  generated  by  the  input 
magnetic  fields.  The  magnetic  field  components  are  the 
desired  results  therefore  the  system  transfer  function  must 
be  determined  so  that  if  can  be  removed  from  the  output 
voltage.  In  the  time  domain  the  relationship  between  the 
output  voltage  and  the  magnetic  field  is  a  convolution 
integral: 

V0i(t)=  Jhi(t-t')  bi(t’)It»  (1) 

where  i=x,y,z.  It  is  more  convenient  to  express  equation  1 
in  the  frequency  domain  by  performing  a  Fourier  transform: 

y=Hi(u>)  Bi**>)  .  (2) 

The  transfer  function  Hi  (u>  )  can  easily  be  determined  in 
terms  of  the  ratio  of  output  to  input. 

The  next  sec+ion  discusses  the  components  of  the 
system  and  attempts  to  develop  a  system  model.  The  last 
section  describes  the  calibration  measurements  and  the 
experimental  determination  of  the  transfer  function. 


PRE-AMPLIFIER 
AMD  FILTER 


PRE-AMP LI  FIE? 
AND  FILTER 


PRE-CONDITIONER 


PRE-CONDITIONER 


PRE-AMPLIFIER 
AND  FILTER 


PRE-CONDITIONER 


Figure  C. 1  System  Component  Configuration 


2.  Svstem  Description 

The  coil  magnetometers  are  composed  of  18  gauge 
copper  wire  with  approximately  5463  turns.  The  individual 
transfer  function  for  the  coil  is: 

VcM»j  KB  (oj)  (3) 

where  Vc  (ujJ  is  the  spectrum  of  the  volrage  directly  out  of 
the  coil  and  K  is  a  gain  factor. 

The  next  component  in  the  system  is  the  preamplifier 
low  pass  filter.  The  preamplifier  essentially  introduces  a 
gain  of  1000  so  the  voitaje  signal  can  be  sent  to  remote 
instrumentation.  The  low  pass  filter  is  a  five  pole 
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Chebyshev  with  a  3dB  cutoff  at  23  Hz.  The  filter  is 
primarily  used  to  remove  53  Hz  sigaaLs  from  the  system  ana 
to  prefilter  the. data  before  it  is  digitalized.  An  addi¬ 
tional  cole  is  introduced  because  of  a  capacitive  coupling 
between  the  second  and  third  stages  of  the  filter.  Figure 
C.2  shows  a  pole-zero  diagram  of  the  preamplifier.  The 
transfer  function  takes  the  fora: 

H(jvj)  (jw+b)  (ju/s-ji)  (jw+c  +  jdi  {-jw+e-jf)  (jui+e+jf)  (4) 

where  K  is  the  gain  of  the  preamplifier  and  the  coeficients 
a  -  f  are  taken  from  the  pole-zero  diagram. 

The  preconditioner  has  additional  gain  of  approxi¬ 
mately  1000  and  has  zener  diodes  to  limit  the  input  voltage 
to  the  cC3  board.  TypicalLy  the  cutoff  is  around  7  to  8 
volt  s. 

The  magnitude  and  phase  response  of  the  preamplifi¬ 
er-filter  and  preconditioner  are  shown  in  Figures  c.3 
through  C.8.  The  input  test  signal  generating  these  curves 
had  fiat  magnitude  and  phase  character  as  shown  in  Figure 
C.9.  The  figures  then  represent  the  actual  transfer  func¬ 
tions  of  the  preamplifier- f liter  preconditioner  system  for 
the  x-y-z  channels  of  land  and  ocean  sensor  systems.  The 
figures  show  that  the  channels  are  weLl  marched  in  magnirude 
and  phase.  The  phase  response  is  linear  and  virtually  iden¬ 
tical  in  every  case.  This  is  an  important  point  since  the 
phase  part  of  the  transfer  function  is  not  removed  from  the 
output  voltage  signal.  However,  since  the  data  analysis  is 
in  terms  of  power  spectral  densities  and  coherences  this 
fact  has  r.o  bearing  on  the  results. 

The  total  transfer  function  of  the  system  is  the 
product  of  equations  3  and  4.  Based  on  these  formulas  and 
the  figures,  the  resulting  transfer  function  should  be 
fairly  linear  at  the  low  frequency  end. 
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Figure  C.3  Eesponse  of  the  Land  X  Systea 


X*  15.781 

TRAMS 


3.  Determination  of  the  System  Transfer  Function 

The  magnitude  of  tie  system  transfer  function  was 
determined  experimentally  through  the  use  of  a  Helmholtz 
coil,  a  source  of  sinusoidal  current  and  a  spectrum 
analyzer.  The  experimental  set  up  is  illustrated  in  Figure 
C.  ID.  A  or.e  turn  Helmholtz  coil  with  a  0.61  meter  radius  and 
a  0.61  meter  separation  between  coils  is  used  to  provide  a 
calibrated  magnetic  field  to  the  sensing  ceil.  The  sensing 
coil  is  placed  concentric  to  the  Helmholtz  coil  and  is 
centered  vertically  to  ensure  a  uniform  magnetic  field  is 
created  about  the  sensing  coil.  The  magnetic  field  in  this 
region  of  the  Helmholtz  coil  is  directed  along  its  axis  and 
is  given  by  [7]: 

Eh=^NI/5;Tr  (5) 


to  terms  on  the  order  of  (L/E)  where  L  is  the  linear  dimen¬ 
sion  of  the  region,  N  is  the  number  of  turns,  .  3  is  the 
Helmholtz  ccrl  radius  and  I  is  the  driving  current  of  the 
coil.  The  driving  current  is  produced  by  a  Wavecek  signal 
generator  model  number  1'42.  The  current  is  measured  by 
reading  the  voltage  across  a  995  ohm  resistor  in  series  with 
the  Helmholtz  coil.  This  voltage  is  measured  by  a  spectrum 
analyzer  (HP-3582A)  to  avoid  any  spurious  signals  generate 
by  harmonics.  Once  the  driving  current  is  known,  then  the 
calibrating  magnetic  field  can  be  determined.  The  voltage 
output  of  the  sensor  system  is  then  measured  by  a  spectrum 
analyzer.  By  repeating  the  measurements  at  a  number  of 
different  frequencies  and  applied  fields,  the  system 
transfer  function  can  be  measured. 

The  actual  calibration  was  conducted  at  the  La  Mesa 
village  site  using  the  setup  illustrated  ir.  Figure  C.10. 
The  Helmholtz  ceil,  sensing  coil  and  preamplifier  are 
located  approximately  40  maters  from  the  shack  which  housed 


Figure  C.  10  System  Calibration  Experiment 


The  remainder  of  the  experiment.  With  the  Helmholtz  coil 
disconnected,  the  background  magnetic  noise  of  this  location 
was  repeatedly  examined  to  ensure  accurate  readings.  Figure 
C. 1 1  shows  a  typical  background  noise  plot  from  the  spectrum 
analyzer  directly  (therefore  the  transfer  function  has  not 
been  removed).  The  distinct  features  are  Schumann  reso¬ 
nances,  otherwise  the  background  was  low  enough  and  stable 
enough  to  adequately  calibrate  the  system.  Only  at  the 
Schumann  frequencies  was  the  background  subtracted  from  the 
measured  output  voltage.  Salibration  runs  were  conducted 
from  0.05  Hz  to  20.0  Hz  with  applied  calibrated  magnetic 
field  of  0.02,  0.2  and  1.0  nanoteslas  (gamma).  Figures  C.12 
through  C. 16  show  the  resulting  calibration  curves  of  the 
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land  and  ocean  magnetometer  systems.  lie  true  transfer  func¬ 
tion  is  represented  by  the  0.02  nl  field  of  the  Helmholtz 
ceil.  The  larger  applied  field  snw  saturation  by  the 
system  caused  by  the  protective  zener  diodes  in  the  precon¬ 
ditioner.  At  the  low  frequency  end,  all  the  systems  are 
closely  matched.  At  1  Hz  an  output  voltage  of  1  v  indicates 


Vo  4(mV) 
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::  ^ 
1 


} 

-  -I 


Figure  C.11  Background  Geomagnetic  Soise  f (Hz>  35/20/82 
151 3  hrs. 

a  measured  magnetic  field  of  0.4  nl.  A  1  nl  field  is 
saturated  above  2  Hz. 

The  computation  of  the  actual  equations  representing 
the  magnitude  of  the  transfer  functions  presented  a  distinct 
problem  because  fitting  the  data  to  tie  theoretical  transfer 
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(volts/nT) 


Figure  C.  12  Land  System  X-Coil  Calibration 

function  cf  the  previous  section  is  difficult.  Ideally  this 
is  desirable  since  the  phase  character  cf  the  systi. 
transfer  function  can  then  be  inferred.  However,  since  the 
priaary  area  of  interest  is  in  the  lower  frequency  region, 
the  transfer  function  was  modeled  in  terms  of  linear 
segments.  The  equations  derived  in  this  manner  have  a 
maximum  error  of  1%  in  the  range  of  3  to  10  Hz  and  5%  in  the 
10  to  20  Hz  range. 

In  the  actual  computation  each  curve  was  broken  into 
convenient  segments  that  could  be  represented  by  a  straight 
line.  The  range  of  the  segnents  were  0  to  3  Hz,  3  to  5  Hz,  5 
to  7.5  Hz,  7.5  to  10  Hz,  10  to  15  3z  and  15  to  20  Hz.  For 
each  segment  the  slope  was  first  computed  by: 


£  | (volts/nT ) 


Figure  C. 13  Land  System  Y-Coil  Calibration 

M=  (  Y2-  Y 1 )  /  ( X2 *  XI ) .  (6) 

Ones  the  slope  of  the  line  vas  found  the  intersection  with 
the  Y  axis  (b)  was  found  from: 

Y=ti  X+b .  (7) 

These  computations  were  performed  for  each  segment  to  yield 
a  series  of  five  equations  to  model  each  sensing  coil  and 
amplification  system.  The  individual  transfer  functions  can 

be  found  in  the  computer  program,  see  Table  I. 
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Figure  C. 1 5  Ocean  System  X-Coil  Calibration 
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Transfsr  Function  3  omputa  r  Program 


SEA: 

1  IF 


IF j FRQ. 
XXO  (L)  = 
YYO  (L)  = 
3  0  TO  8 
IF  (FRQ. 
XXO  (L)  = 
YYO  JL)  = 
30  TO  8 
IF  (FRQ. 


LE.15.)  GO  TO  2 

XXO  (L)  /  (-8  1 .36*17.  066*  FRQ- 3.  5  389*  FRQ* FR  2) 
YYO  (L)  /  (-63.435  +  11  .7106*FRQ- 0. 23279*FEQ*F] 


LE.10.  )  GO  TO  3 

XXO  (L)  /  (-18.486  +  4.  9029  *FRQ) 

YYO  (L)  /  (-6  3.435  +  1  1 .710 6* FRQ- 0 . 23279* FRQ* FRQ) 

LE.7.5) GO  TO  4 

XXO  (L)  /  (3.  947*FRQ-  9.35  81 

YYO(L)  /  (4.  295*FRQ- 11 .837) 


IF (FRQ. 
XXO  (L)  = 
YYO  (L)  = 
3  0  TO  8 


LE.5.)  GO  TO  5 

XXO  (Li  /  (2.  8  105* FRQ -0.5  6  57) 

YYO  (L)  /  (2.  8  00*FRQ  -  .50) 


5  XXO 

(L)  =XXO(] 

-)  /  ( 

YYO 

(l)  =yyo(i 

-)  /  ( 

7  17*FRQ» 


LAND: 

1  IF  (FRQ 
XXX  (L) 
YYYjL) 
30  TO 

2  IF  (FRQ 
XXX  (L) 
YYYjL) 
30  To 

3  IF  (FRQ 
XXX  (L) 
YYYjL) 
3  0  TO 


4  I F (FRQ 
XXX  (L) 
YYYjL) 
3  0  TO 

5  IF  (FRQ 
XXX  (L) 

7  YYYjL) 
3  0  TO 

6  XXX  (L) 


.  LE.15. )  GO  TO  2 
=  XXX(L)  /  (10  5.5-3. 14  *FRQ) 

=  YYY  (L)  /  (181.32-7.  388*?RQ) 

8 

i.  LE.10.  )  GO  TO  3 

=XX X  (L)  /  (5.  958*FRQ- 30.9  7) 

=  YYY  ( L)  /  (7.  166*FRQ-  39.9  9) 

8 

.  LE  .7.5)  GO  TO  4 
=  XXX(L)  /  (3.  492*FRQ-  6.  31) 
=YYY(L)  /  (4.  252*FRQ-  10.3  5) 

8 

.  LE.5  .)  GO  TO  5 
=  XXX(L)  /  (2.  6311*FRQ+0.14667) 
=  YYY  (L)  /(3.012*FRQ-1.55) 

8 

i.LE.3.)  GO  TO  6 
•  XXX  (L)  /  (2.  6311*FRQ  +  0.1  4667) 
=  YY Y  (L)  /  (2.7  02*FRQ) 

=XX X  (L)  /  (2. 72*FRQ) 
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APPENDIX  D 


A.  PDLSE  CODE  MODULATION  SISTEM 

The  PCM  circuit  compensate  are  minted  on  a  3  inch  by  15 
inch  printed  circuit  board  with  a  15  pin  double  sided  edge 
connector  at  each  end.  The  pia  assignments  for  the  signal 
connector  and  the  control  connector  ace  listed  ia  tables  D.1 
and  D. 2,  respectively.  The  circuit  is  diagramed  in  figures 
A. 4  and  A. 5. 


TABLE  II 

Signal  Connector  Pin  Ass:  ^nment 


PIN 

"I" 

B 

C 

D 

E 

F 

H 

J 

K 

L 

S 

N 

4» 

H 

S 

1 

3 

15 


ASSIGNMENT 
cITan ne I  T  i n o 1 1 
channel  2  inoat 
channel  3  input 
channel  4  inoat 
channel  8  input 
channel  7  inoat 
channel  6  input 
channel  5  input 
channel  15  inbut 
channel  14  inbit 
channel  13  input 
channel  12  inoit 


channel  11  input 
channel  10  inout 
channel  9  input 
-157  (max.) 
common  (signal) 

♦  15v  (max.) 


analog) 
analog) 
a  nalogi 
analog) 
a  nalogi 
(analog) 
(a  nalogi 
analog) 
a  nalogi 
analog) 
a  nalogi 
analog) 
a  nalogi 
analog) 
a  nalogi 


The  PCM  system  incorporates  a  crystal  oscillator  and 
associated  CMOS  integrated  circuitry  to  develop  the  clocking 
pulses,  a  16  channel  CMOS  analog  multiplexer,  a  12  bit  CMOS 
analog-to-digital  converter  and  associated  circuitry  to 
provide  the  pulse  coding. 
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1 


2 

3 

4 

5 

6 

7 

8 
9 

10 
1  1 
12 

13 

14 

15 
A 


D 

c 

D 

K 

L 

a 

N 

P 

R 

S 


TABLE  III 

Control  Connector  Pin  Assignment 


SIGNMENT 

TU¥T"ITI  (CHSL  1 4| 


WO 

EXTDW2  IN 
EXTDW2  SEL 
SAMP.  RATE  SEL  IN 
DCL 

POWER  GRD. 

♦  157 
NRZ  -  L 


Bi0-L 
EXT  ffi 
Cl  t 
C2 
C  4 
C  8 

TX  FILTER  INPUT 
S 13  GND 

TX  FILTER  OUTPUT 

SIS  GND 

-157 

+127  INPUT 
SAMPLE  RATE  123 
SAMPLE  RATE  64 
SAMPLE  RATE  32 
SAMPLE  RATE  15 
SAMPLE  RATE  8 


The  crystal  clock  oscillator  operates  at  a  frequency  of 
24.576  KKz  producing  a  square  wave  output  with  a  loss  bit 
rate  of  one  bit  in  one  million.  This  pulse  train  with  pulse 
duration  cf  20.35  microseconds  triggers  the  clock  input  of  a 
binary  up  counter  Z17(1/2i  and  is  aLso  routed  directly  to 
pin  M  of  the  control  connector.  Pin  2  of  Z17  is  the  counter 
enable  input  which  is  connected  to  7+  thereby  constantly 
enabling  the  counter.  Pin  7  is  the  reset  input  which  is 
grounded.  The  counter  is  composed  of  four  toggle  flipflops 
that  are  incremented  on  the  positive  joing  edge  of  the  clock 
pulses.  The  flip-flop  outputs  Ql  through  Q4  are  routed  to 
control  connector  pins  N,P,R,  and  S,  respectively.  The  pulse 
repetition  rate  at  pins  M,N,P,  3,  aid  S  represent  sample 
rates  of  128,  64,  32,  16,  and  8.  Selection  of  the  sample 
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rats  is  dependent  upon  the  Syguist  sampling  criterion.  If  f 
is  the  highest  frequency  of  Interest  in  the  analog  data 
being  sampled,  than  the  sample  rata  should  be  equal  to  or 
greater  than  2f  for  unbiasad  resolution.  Currently  a  sample 
rata  of  64  par  second  is  employed  to  ensure  the  Ryquist 
sampling  criterion  is  mat  ovar  the  frequency  range  of 
interest  (0-20  Hz).  A  connection  between  the  pin  that 
represents  the  selected  sampling  rata  and  pin  5  of  the 
control  pin  connector  will  establish  that  pulse  repetition 
rata  as  the  system  clock. 

The  clock  triggers  Z9,  a  four  stage  binary  up  countar. 
Z8  inputs  consist  of  a  single  clock,  carry-in  (clock 
enable),  binary/decade,  up/do«n,  preset  enable,  and  four 
individual  jam  signals.  Four  seperate  buffered  2  signals  and 
a  carry  out  signal  are  provided  as  outputs.  V*  input  on 
pins  9  ar.d  10  configure  the  circuit  as  an  up  binary  counter. 
The  jam  signal  is  selected  to  develop  the  data  word  length. 
The  counter  starts  at  a  count  of  4  as  configured  and  counts 
to  15  (12  counts).  CCO  goes  high  setting  Z12.  The  preset 
enable  signal  on  Z8  receives  the  high  set  output  from  Z12 
(STIR)  and  the  Z 12  counter  is  asyncaronously  preset  to  the 
jam  signal  count.  The  preset  count  is  gated  through  Z9  and 
Z10  and  the  convert  start  signal  to  Z3,  the  16  channel 
analog-to-digital  converter.  The  carry  out  signal  (CCO) 
from  the  Z8  counter  acts  as  the  clock  input  to  Z17(l/2) 
another  four  stage  binary  up  counter.  Its  outputs,  Q1 
through  Q4  represent  the  multiplex  channel  enable  signals  to 
Z2,  the  16  channel  analog  multiplexer.  Each  channel  is 
enabled  for  12  clock  puLses.  The  analog  multiplexer, 
depending  on  the  state  of  tae  multiplex  enable  signal, 
switches  a  common,  output  to  one  of  sixteen  inputs.  Its 
output  is  amplified,  conditioned  and  passed  to  the  19  Volt 
input  of  Z3,  the  analog-to-digital  converter.  The  parallel 


digital  output  of  Z3  is  passed  to  Z4  and  Z5  for  paraliel-to- 
serial  conversion.  The  analog  iata  word  (ADW)  is  gazed  to 
NANO  gate  Z1 1  along  with  Wo,  D CL,  aid  EXTWD2.  Wo  and  Wo  are 
developed  in  NOR  gate  Z9(1/2).  Wo  is  high  when  Cl  through  C8 
are  zero  and  low  otherwisa.  Then  using  Z17*s  counter  zero 
count  as  intial  sequencing  point,  the  frame  sync  word  is 
passed  through  Z  13  pin  13  followed  by  channel  1  through 
channel  15  analog  data  words.  Z12  and  Zlt*  develop  the 
signals  to  be  biphase  logic  or  non-return  to  zero  (NR Z) 
logic.  This  system  feeds  the  biphase  logic  BI^-L  through  a 
low  pass  filter  and  a  signal  conditioning  circuit  where  the 
modulation  index  is  set.  tn  the  sea  system  the  offset  and 
gain  potentiometers  are  adjusted  so  the  signal  conforms  to 
TTL  levels.  This  final  output  is  the  iata  input  to  the  optic 
tran  smitzer. 
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A.  OPTICAL  FIBERS 
1 .  Fiber  types 

Although  many  types  of  fibers  have  beer,  cried  in  the 
recent  past,  present  technology  favors  three  classes  of 
fibers  for  optical  communications:  single-mode  step-index, 
multi-mode  step-index,  and  multi-mode  graded- index  (see 
Figure  E.  1) . 

a.  Single-Mode  Step-Index 

Single-mode  fibers  are  tyolified  by  a  very  thin 
core  (1.5  to  8  micrometers)  with  uniform  index  of  refraction 
surrounded  by  a  cladding  of  considerably  greater  thickness 
whose  index  of  refraction  differs  by  a  very  small  amount 
from  that  of  the  core.  The  "single  mode"  feature  can  be 
attributed  to  the  abrupt  interface  between  core  and  cladding 
as  well  as  the  small  change  ia  the  refractive  index.  The 
cladding  material  is  lossy  so  that  any  modes  it  accepts  will 
guickly  attenuate  and  not  couple  back  into  the  core.  Since 
the  difference  in  refractive  index  from  core  to  cladding  is 
small,  the  numerical  aperture  is  small  also  on  the  order  of 
0.14.  The  small  aperture  makes  coupLing  to  source/detector 
difficult  [8,  9]. 

b.  Multi-Mode  Step-Index 

The  core  diameter  of  multi-mode  step-index 
fibers  are  larger  (25-150  micrometers).  The  refractive 
index  is  uniform  across  the  core  diameter  [9],  The  cladding 
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refractive  index  is  slightly  lower  than  that  of  the  core. 
The  numerical  aperture  is  small  but  the  core  size  is  larger 
than  that  of  the  single  mode  fiber,  aLleviating  somewhat  the 
coupling  problem.  The  loss  characteristics  of  these  fibers 
can  be  extremely  low  and  in  general  tan  be  linked  directly 
to  the  lower  numerical  apertures  and  manufacturing  tech¬ 
niques  employed. 

c.  Multi-Mode  Sraiad-lndex 

Graded  index  fibers  appear  to  be  the  optimum 
fiber  (of  those  currently  available!  for  long  line  applica¬ 
tions  (8].  The  fibers  have  core  diameter  ranging  from  20  to 
150  micrometers  with  an  index  of  refraction  that  varies  as  a 
function  of  the  distance  from  the  core  center  as  shown  in 
Figure  E.2.  The  grade  in  the  index  removes  the  abrupt  change 
that  is  typical  of  the  step- index  fioar.  The  net  effect  is 
that  all  modes  cross  the  optical  axis  at  about  the  same 
places  and  at  about  the  same  time  [ID]. 

2.  Attenuation  in  Fiber  Wav equides 
a.  Losses  in  General 

The  principal  causes  for  attenuation  losses 

are: 

1)  Material  absorption:  This  loss  can  be  mostly 
attributed  to  Hydroxyl  and  Transition  metal 
ions  in  glass. 

2)  Scattering:  This  loss  Is  primarily  caused  by 
geometric  irregularities  at  the  core-cladding 
in  terf  ace. 

3)  Radiation  losses:  This  Loss  is  a  result  of 
bending  the  fiber  cable  a r  small  radii  of 
curvat ure  [11]. 


Scattering  Losses 


b. 


Ml  transparent  materials  scatter 
due  to  frozen  in  thermal  fluctuations  of  constituent  atoms. 
This  intrinsic  scattering  represents  the  fundamental  contri¬ 
bution  to  attenuation  in  optical  waveguides.  The  glasses 
utilized  in  fiber  waveguides  generally  have  more  than  one 
oxide  present  that  cause  another  form  of  scattering.  This  is 
due  to  concentration  fluctuations  in  the  constituent  oxides. 
For  high  silica  glasses,  concentration  fluctuations  typi¬ 
cally  account  for  only  approximately  25*  of  the  observed 
scattering  loss. 

In  addition  to  these  two  intrinsic  scattering 
losses,  induced  scattering  through  non-linear  effects  can 
exist  such  as  stimulated  Raman  and  Brillouin  scattering. 
Because  of  the  snail  core  size,  the  confined  guidance  and 
the  long  interaction  length,  relatively  low  absolute  power 
levels  are  required  to  observe  such  effects.  Actually  2000 
watts  has  been  injected  into  a  2  Km  length  of  guide  having  a 
75  micrometer  core  diameter  with  no  non-linear  attenuation 
obse  rvea  (11]. 


c.  Radiation  Losses 


These  losses  are  associated  with  the  waveguide 
structure.  The  cladding  has  in  practice  been,  a  thickness  of 
a  few  tens  of  micrometers.  If  the  jacket  is  lossy  (to 
minimize  crosstalk),  some  fraction  of  the  mode  field  can 
reach  this  region  and  be  attenuated.  Additional  radiation 
losses  are  experienced  when  the  optical  wave  guide  is 
curved.  The  presence  of  outside  perturbations  on  a  fiber 
cause  light  loss  through  the  cladding.  In  step-index  fibers 
this  mechanism  causes  a  change  in  the  angle  of  incidence  of 
the  light  rays  until  a  critical  loss  angle  is  reached.  In 
graded-ir.dex  fibers  loss  occurs  from  mode  mixing  and  loss  of 


the  higher  order  modes  in  a  simular  fashion.  These  losses 
can  become  appreciable  in  fiber  cables  for  underwater  appli¬ 
cation.  Losses  as  great  as  U3Q  dB/Km  have  been  reported 

r  i2i. 

d.  Loss  Minimization 

Fortunately  it  is  known  than  fibars  can  be 
protected  from  microbendiag  and  encased  in  constructions 
capable  of  high  pressure  [13].  The  size  of  these  losses  is 
typically  about  1  dB/ka.  Minimization  of  microbeading  losses 
is  achieved  by  encasing  the  fiber  in  alternately  soft  and 
hard  jackets.  A  transverse  force  applied  to  the  stiff  outer 
jacket  will  be  resisted  and  spread  over  a  large  area.  The 
soft  inner  jacket  will  comoLy  with  the  force,  attenuating  it 
and  spreading  it  ever  a  still  larger  area.  The  total  force 
will  then  be  able  to  transmit  to  the  fiber  which  is  itself 
stiff.  For  example  FTT's  T— 1211  fiber  is  buffered  with 
8TV184,  a  silicon  rubber,  during  the  drawing  process.  This 
has  the  advantage  of  protecting  the  fiber  from  contaminants 
early  on,  particularly  QH  ions.  A  second  buffer,  Hytrel 
7245,  a  hard  polyester  elastomer.  Is  then  applied.  This 
combination  makes  the  fiber  a  good  choice  for  underwater 
cabling.  Hytrel  also  gives  the  fibers  hundreds  of  hours  of 
immunity  from  moisture  [13]. 

e.  Losses  due  to  Hydrostatic  Pressure 

Hydrostatic  pressures  can  also  affect  the  atten¬ 
uation  in  optical  fibers.  Silicon  rubber  clad  silica  (PCS) 
is  compressed  in  sea  water.  At  2000  psi,  the  numerical  aper¬ 
ture  (NA)  for  this  type  fiber  approaches  zero  and  no  light 
is  transmitted  due  to  compression  of  the  silicone  and 
increase  of  the  s  licoa?  index  of  refraction  [14], 
Therefore,  glass  claf  ibers  must  be  used  in  underwater 
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3*  Static  Fatigue  and  5  trena th 

a.  Since  ihs  fiber  optic  utilized  in  a  data  link 
for  magnetic  data  collection  need  not  be  a  supporting 
element  of  the  system  the  tensile  strength  is  not  a  signifi¬ 
cant  factor. 

While  the  theoretical  strength  of  optical  fiber  is 
more  than  600  kg/sq  mm  *15]  the  actual  strength  is  related 
to  minimum  strength  points  due  to  flaws  and  cracks  on  the 
fiber  surface.  For  a  breaking  strain  of  1  percent,  the 
minimum  flaw  size  must  ne  lass  than  1  micrometer,  which  is 
difficult  to  assure  in  production  rune  "16].  Further  more, 
the  growth  of  flaws  is  directly  related  to  temperature, 
applied  strain,  relative  humidity  and  even  pH.  \  pH  of  7-8 
causes  the  most  severe  flaw  growth  *17].  This  value  approxi¬ 
mates  the  pH  of  the  ocean. 

4.  Cable  Con  figuration 

Cabling  cf  opticaL  fibers  iae  the  objective  of 
providing  useable  strength  at  minimun  weight  while  ensuring 
environmental  protection.  Extensive  tests  upon  several 
different  configurations  of  mini-cabLes  for  underwater  use 
with  and  without  conductors  has  been  conducted  by  NOSC, 
Hawaii  [18].  The  following  presents  some  of  these  cables 
whose  characteristics  may  make  them  suitable  for  underwater 
magnetic  data  collection  operations. 

Figure  5. 2  represents  three  possible  cable  arrange¬ 
ments,  all  similar,  except  for  the  composition  of  the  load 
bearing  structure.  This  may  be  ail  Kevlar,  a  hybrid  mix  of 
Kevlar  and  S-Glass,  or  all  3-Glass  in  an  epoxy  matrix. 

a.  Lightweight  Single  Fiber  Sable 

1)  Cable  5358  (all  Kevlar):  This  cable 
survived  500,000  flexure  cycles  at  23  percent  loading.  The 
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Figure  E. 2  Lightweight  SingLa  Fiber  Cable 
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breaking  strength  of  this  cable:  is  223  pounds  with  a  tensile 
strength  of  130,  2  00  psi. 

2)  Cable  5033  (Hybriil  :  This  cable  survived 
2,370, 0G0  flexure  cycles  without  failure  at  20  percent 
loading.  Its  breaking  strength  is  227  pounds  and  tensil® 
strength  is  13  0,  0  00  psi. 

3)  Cable  3033  (all  5-3lass)  :  This  is  the 
strongest  combination  virh  a  breaking  strength  of  346  pounds 
and  a  tensile  strength  of  275,000  psL.  This  caole  survived 
750,000  flexure  cycles  at  20  percent  loading  without 
breaking  before  the  test  was  stopped. 

Those  cables  tasted  by  113 3 C  Hawaii  provide  the 
extension  of  fiber  optic  technology  to  the  underwater  envi¬ 
ronment  and  those  mentioned  above  allow  for  utilization  in 
underwater  magnetic  data  collection. 

B.  LIGHT  SOURCES  AH D  DETECTORS 

Light  sources  for  fiber  optic  systems  require  certain 
characteristics  including  long-life-time  in  use,  high  effi¬ 
ciency,  reasonably  low  cost,  sufficient  optical  power 
output,  capability  for  various  types  of  modulation  and 
physical  compatibility  with  fiber  ends.  There  are  three 
basic  configurations:  light  emitting  diode,  solid  state 
laser,  and  semiconductor  injection  laser. 

1.  Light  Sources 

a.  Light  Emitting  Diode  (LSD) 

LED's  generate  light  wnen  carriers  injected 
across  a  p-n  junction  radiatively  recombine,  emmitting  radi¬ 
ation  over  a  solid  angle  of  2ir  steradians.  The  actual 
output  power  emitted  into  a  unit  solid  angle  per  unit  emis¬ 
sion  area  is  low  for  LED's  even  though  the  total  power 
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output  is  reasonably  higi.  Thais  are  two  types  of  LED 
sources  available,  edge  emitters  and  domed  surface  emitters. 
Domed  emitters  are  rather  complex  to  manufacture  resulting 
in  edge  emitter  LED's  being  mere  commonly  utilized.  Noise 
introduced  by  LED's  into  tne  optical  fibers  is  relatively 
small  in  most  cases  and  is  proportional  to  the  spectral 
width  and  the  number  of  transmission  modes  propagated  in  the 
fiber.  Typical  power  verses  drive  response  curves  for  LED's 
are  usually  linear  (see  Figure  E.3>.  dost  available  LED 
sources  have  peak  wavelengths  of  32)  to  850  nm  and  power 
outputs  of  approximately  1  nv.  LED's  are  constructed  to  last 
in  excess  of  1  million  hours  [19,  S]. 


OO  iso  200  256  300  390 
DIRECT  CURRENT  InA) 


Laser*  DH,  stripe  width 
quantum  eff.  3055 

t  Homo  junction, 
Dome  emitter, 
.4nun  dia 


Figure  E.3  Typical  Power  vs  Drive  Response  Curves 


b.  Lasers 

Double-ketero  junction  DH  Laser  structures  car.  be 
made  with  relatively  long  lifetimes  when  the  lattice  parame¬ 
ters  of  the  electrical  composition  in  the  various  layers  and 
boundaries  are  carefully  Hatched.  In  the  DH  laser  with 
stripe  geomerry,  the  lasing  region  occurs  only  in  the  trans¬ 
verse  direction  under  the  width  of  the  contact  stripe. 

The  basic  semico nduc tor  material  for  both  LSD's 
and  injection  lasers  is  the  GaAs  combination.  Doping 
elements  like  Al ,  As,  In,  5e,  and  Sb,  are  used  to  shift  the 
emission  wavelength.  The  threshold  current  for  typical  DH 
lasers  range  from  0.7  to  14  kA/sq  cm  and  the  quantum  effi¬ 
ciency  from  10  to  50*?  compared  with  3%  for  an  LED  resulting 
in  a  large  incremental  increase  in  output  power  for  small 
increases  in  current  beyond  -.ae  current  threshold  (Figure 
E.4»  . 

c.  Comparison 

Virtually  all  applications  of  fiber  optic  tech¬ 
nology  for  local  data  communications  to  date  have  chosen 
LED's  for  use  as  the  sourcn.  Bandwidth  distance  rates  have 
been  within  the  capabilities  of  this  device  and  its  cheaper 
cost,  ease  of  use  and  operational  expediency  have  made  it 
competitive  with  lasers. 

2.  Derectors 

In  fiber  optic  systems,  the  most  commonly  used 
receivers  utilize  photodiodes  (either  PIN  or  avalanche 
types)  ro  convert  incident  Light  into  electrical  energy. 


a 


PIN  Photodiode  (PIN)^ 


The  PIN  diode  consists  of  a  relatively  large 
(very  lightly  doped)  region  sandwiched  between  d  and  n  doped 
semiconducting  reaions.  Pnotons  absorbed  in  this  region 
create  electron-hole  pairs  that  are  then  separated  by  an 
electric  field,  thus  generating  an  electric  current  in  the 
load  circuit.  The  quantum  efficiency  of  the  photodiode  is  a 
function  of  wave  length  and  temperature  [8]* 

b.  Avalanche  Photodiodes  (APD) 

The  APD  is  designed  for  applications  requiring 
greater  sensitivity.  Because  of  a  strong  electric  field 
arising  within  it  as  a  result  of  external  biasing,  the  APD 
exhibits  an  internal  gain  nechanisa.  They  require  a  consid¬ 
erably  higher  bias  voltage  than  PIN  diodes  (bias  voltages  on 
the  order  of  300  V  are  not  urcomr.cn)  *3]. 

t 

c.  Comparision  of  PIN  and  API 

As  a  result  of  the  Ugh  reverse  bias  require¬ 
ments  of  the  APD  and  the  power  constraints  of  the  magnetic 
data  collection  system  che  PIN  diode  is  the  obvious  choice 
for  the  system's  detector. 

C.  TR  A  NS  HITT  EH/ RECEIVES 

There  is  a  wide  variety  cf  c oumercially  available  trans- 
mitter/receivsr  sets  wit a  tne  full  range  of  scarce/ detectors 
previously  mentioned.  These  preclude  the  necessity  of 
designing  the  associated  circuitry.  The  Burr-Brown  3713T/R 
is  such  a  set. 


D.  SPLICING  AND  CONNECTING 


A  problem  in  many  applications  is  in  joining  one  length 
of  fiber  or  cable  to  another,  either  ny  splicing  (relatively 
penanent)  or  using  collectors  (intended  for  repeated 
disconnection  and  reconnection}. 


1 .  Splicing 


Improvements  in  t anuf actar in g  techniques  that 
continue  to  lower  fiber  attenuation  aave  increased  the  need 
for  lew-loss  field  splicing  methods  tiat  will  permit  greater 
spacing  between  repeaters.  The  coupling  loss  of  a  splice 
depends  on  two  types  of  factors:  those  that  are  intrinsic  to 
the  fibers  being  spliced  and  extrinsic  ones  which  depend  on 
the  splicing  technique.  Intrinsic  factors  rncluie  core  dia¬ 
meter,  numerical  aperture,  and  refractive  index  profile. 
These  can  be  expected  to  vary  randomly  because  they  are 
influenced  by  the  manufacturing  and  quality  control 
processes.  Their  effect  on  splice  loss  can  be  treated 
statistically.  The  factors  introduced  oy  the  splicing  method 
include  lateral  m  isalignmen t ,  end  separation,  and  end  prepa¬ 
ration  quality.  Of  these,  lateral  misalignment  is  the  most 
critical. 

The  ideal  splicing  metned  is  a  field  usable  tech¬ 
nique  which  minimizes  these  fiber  extrinsic  factors. 
Splicing  methods  now  avaiiible  depend  on  either  fusion  or 
adhesives  f 2 0  1  - 

a.  Fusion  Method 


The  fusion  or  "hot  splice"  method  uses  an  elec¬ 
tric  arc  or  other  heat  source  to  fuse  two  fiber  ends 
together  after  they  have  been  accurately  aligned. 
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b.  Adhesive  Method 


In  the  adhesive  Dr  mechanical  approach,  a 
mechanical  device  is  used  to  align  the  fibers  through  refer¬ 
ence  to  their  outer  diameters.  The  fibers  Then  are  bonded 
into  this  position  permanently,  usualLy  with  an  epoxy 

adhesive. 

2.  Connectin  a 

The  same  intrinsic  and  extrinsic  factors  apply  to 
the  design  of  fiber  optic  connectors  with  the  added  require¬ 
ment  that  the  connector  system  permit  numerous  disconnec¬ 
tions  and  rematings  without  an  appreciable  increase  in  the 
coupling  loss.  There  are  two  basic  types  of  connectors: 
those  which  take  the  alignment  reference  from  the  outside 
diameters  of  the  fibers  themselves,  aid  those  which  encapsu¬ 
late  the  fibers  then  reference  the  outer  diameters  of  the 
encapsulating  members. 

a.  Capillary  Type 

The  most  obvious  design  ossed  or.  alignment  of 
the  fibers  themselves  is  the  mating  of  two  fibers  inside  a 
precision  capillary  with  diameter  only  slightly  greater  than 
that  of  the  fibers.  A  variation  of  the  capillary  approach  is 
the  quasi-capillary  formed  oy  three  precision  rods. 

b.  Other  Types 

The  other  commercial  terminations  fail  into  the 
second  category.  Some  use  a  molded  taper  or  cylindrical 
ferrule  to  align  one  fiber  with  another.  Most  of  the  avail¬ 
able  connectors  which  use  mechanical  alignment  of  machined 
members  are  "dry”  (they  do  not  use  an  index  matching  fluid), 
and  many  are  based  on  the  5MA-typa  connector  developed  for 
radio  frequency  transmission  over  copper  coaxial  cable  [  20  ]. 
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E.  OPTICAL  DATA  LINK  EN3INE  ERIS  G 


The  following  paragraphs  outline  procedures  for  ergi- 
neecing  a  fiber  optic  data  link  utilizing  the  3urr-3rown 
3713T  transmitter  and  3713R  receiver  [21]. 

1 .  Determining  the  Table  Typt 

Initially,  the  primary  -factor  which  determines 
subsequent  design  of  a  fiber  optic  system  is  tie  amount  of 
attenuation  expected.  Fibers  have  been  developed  which 
exhibit  attenuation  of  less  than  1  dB/km.  Commercial, 
fibers  are  available  with  attenuation  of  less  than  6  dE/km 
at  3.85  micrometer  [22].  lie  ruggedized  graded-index  optical 
cables  that  were  tested  by  NOSC ,  Hawaii  had  attenuation  of 
approximately  7  dB/km  [23].  In  most  applications  the  cost  of 
the  fiber  is  inversely  proportional  to  its  attenuation.  For 
this  reason,  selection  of  a  fiber  is  more  directly  related 
tc  dollars  than  dB  [22].  There  are  also  environmental 
concerns.  The  mechanical  stresses,  environmental  temperature 
extremes,  and  exposure  to  narsh  chemicals  must  be  considered 
when  choosing  a  particular  cable.  Such  things  as  docketing 
material  and  the  type  of  added  strength  members  (see  Figure 
E.  2)  determine  a  cable's  mechanical  aid  environmental 

capabilities. 

There  are  three  basic  types  of  optical  fibers  to 
choose  from:  step  index,  single-mode,  and  graded-indax.  Each 
type  specifies  the  profile  or  variation  of  the  fiber's  index 
of  refraction  as  a  function  of  radial  distance  from  the 
fiber's  center  as  previously  diagrammed  in  Figure  E. 1. 
While  step-index  fiber  with  a  core  diameter  of  200  microme¬ 
ters  or  larger  is  recommended  for  use  with  the  3713T  and 
3713R  as  it  is  the  best  compromise  between  core  diameter, 
bandwidth,  coupling  loss,  and  cost,  the  graded  index  type 
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5083  (Hybrid)  cable  previously  described  is  incorporated 
into  the  subject  system  design.  The  outside  ciaaeier 
disparity  is  reduced  by  application  of  several  jackets  of 
heat  shrink  tubing.  For  short  lengths  of  fiber,  the 
increased  coupling  loss  associated  with  the  smaller  diameter 
core  is  tolerable. 

2.  Determining  Launched  Power 

Once  the  mechanical  characteristics  and  the  cab! -a 
type  have  been  chosen,  it  becomes  necessary  to  analyze  lir.k 
performance  with  various  cables.  The  first  analysis  to 
undertake  is  the  power  launched  into  a  particular  cable  from 
the  optic  transmitter. 

Beginning  with  the  LSD  soiree,  four  import?,  vt 
optical  parameters  must  be  known  in  order  to  calculate  accu¬ 
rately  the  device's  interaction  with  the  system.  These  are: 

a.  Total  optical  power 

b.  Peak  wavelength  of  emission 

c.  Physical  size  of  output  port 

d.  Spatial  profile  of  the  emission  pattern 

Total  output  power  (Po)  is  usually  given  in  units  of 
microwatts.  For  convenience  the  optical  power,  as  well  *s 
the  losses  throughout  the  system  and  the  detector  sensi¬ 
tivity,  can  be  expressed  in  dB.  The  total  system  loss  is  the 
sum  of  these  individual  components.  To  convert  tie  LED  total 
power  output  to  dBm,  refer  it  to  a  level  of  ore  milliwatt. 

LED  output  power  (dBm)  =  10  log(Po(uw)/1000l  (1) 

The  peak  wavelength  of  emission  ^  is  necessary  for 
determining  the  loss  of  the  system  fiber  (which  is  a  func¬ 
tion  of  wavelength).  It  is  not  unusual  for  fiber  attenua¬ 
tion  to  double  or  even  triple  over  only  a  103  nanometer 
range  in  wavelength. 


1 1 


5 


areas  of 


A  power  loss  can  occur  if  there  is  a  mismatch  of 
the  emitting  port  and  the  fiber  end.  This  loss 
expressed  in  dB,  is: 

LED-to-fiber  area  los3(dB)  =  201og(D2/D1)  (2) 

where  D1  and  D2  are  the  diameters  of  the  system  cores  and 
the  LED  output  port,  respeot  ively .  If  D1  is  greater  than  D2, 
the  area  loss  is  zero. 

The  spatial  emission  pattern  of  the  LSD's  output 
umst  be  known,  because  a  loss  can  oocur  if  the  high-crder 
mode  rays  of  the  LED's  output  cannot  be  captured  into  the 
numerical  aperture  (NA)  of  the  system  fiber.  This  is  called 
NA  loss.  Most  data  sheets  give  information  on  numerical 
aperture.  If  the  NA  of  the  LSD  is  less  than  that  of  the 
fiber  then  no  NA  loss  occirs  at  that  interface.  If  the 
system  cable  has  the  smaller  NA,  then  the  loss  is: 

NA  loss  (dB)  =  201cg(NA  source/NA  system  fiber)  (3) 

If  there  are  no  fiber  splices  or  taps  in  the  system 
the  remaining  losses  will  os  at  the  fiber  detector  inter- 
face.  These  losses  may  be  calculated  as  described  above. 

Typical  gap  and  misalignment  displacement  in  the 
various  connector  systems  available  for  the  fiber  optic 
active  component  package  are  3.15  mm  and  0.05  mm  respec¬ 
tively.  Angular  misalignment  in  these  connectors  is  very 
small  and  can  be  ignored. 

Beginning  at  the  detector  and  working  backwards,  the 
required  LED  drive  current  can  be  calculated  as  follows. 

Generally  receiver  sensitivity  is  a  function  of 
bandwidth  and  desired  bit-e  rror-r ate  (3ER)  .  The  detector 
data  sheet  specifies  the  minimum  power  to  be  launched  into 
the  detector  input  port  for  a  given  BER.  Therefore,  it  is 
necessary  to  calculate  the  coupling  loss  from  the  system 


fiber  to  the  detector.  This  will  yield  the  optical  power 
required  at  the  detector  end  of  the  fiber.  Ir  is  -;*.en  a 
simple  matter  to  subtract  the  loss  of  the  cable  to  arrive  at 
the  required  launched  power.  Detector  minimum  power  can  be 
converted  to  dBm  by  equation  (1).  Cable  loss  is  calculated 
by  multiplying  cable  length  (km)  by  cable  attenuation 
(dB/km) . 

PL  =  PDM  -  SAID  -  ALD  -  CL  -  FL  -  3dB  (4) 

where, 

PL  =  launched  power 

PDM  =  minimum  detector  power 

SALD  =  detector  SA  loss 

ALD  =  detector  area  lose 

CL  =  cable  loss 

FL  *  Fresnel  loss 

For  design  ease  the  37131/3  specif icatlon  sheet  includes  the 
graph  illustrated  in  Figure  E.  4.  Power  launched  may  be 
determined  simply  by  entering  the  graph  with  the  core 

diameter. 

3.  Determining  Ka xim uja  Cabi§,  Length 

or.ce  the  actual  launched  power  into  a  specific  cable 
(PL)  and  the  attenuation  of  a  fiber  at  a  specific  wavelength 
(L  i  are  known,  the  maximum  cable  length  may  be  obtained. 
First  the  loss  margin  (LM)  is  found  by: 

LM(d3)  =  10  log  { P  L/Pi  n,  mini  (5) 

Pin, min  is  the  minimum  rated  output  power  (with  power 
adjusted  set  for  maximum  power)  usually  in  microwatts. 
Pin, min  is  specified  in  the  transmitter  specification  sheet. 
Then  Xmax,  the  maximum  cable  length  that  will  just  present 
Pin,  min  to  the  transmitter  is  found  by: 
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Figure  2.4  Launched  Power  vs.  Core  Diaaeter 


Xmax  =  LM  (d3) /tyd  Bra)  (6) 

If  this  length  is  too  short,  a  cable  with  a  larger  diamete: 
(larger  PL)  or  less  attenuation  (smaller  L^)  or  a  trans 
nitter  with  a  higher  Pin, rain  must  be  used. 

The  type  5088  cable,  for  exaraple,  has  L^  =  7dB/km  a 
660  nm  and  a  core  diameter  0.  005  inch  or  127  micrometer 
which  corresponds  to  0.9  microwatts.  Pin, rain  =  30  nanowatt 

-9 

worst  case  for  10  BER.  Thus: 

LM  *  101og(0.9  microwatts/30  nanowatts)  (7) 

LM  =  1  0  lo g 3  0  =  14.  77dB  (8) 

Xraax  =  14.  77dB/7x1  0  3dBra  (9) 

Xmax  =  2.11  km  (10) 


This  length  is  longer  than  required  and  is  therefore 
sufficient. 

4.  Determining  Minimum  Cable  Lenqth 

Short  cable  lengths  with  large  diameter  finers  may 
cause  some  receiver  slew  rate  limit  which  will  appear  as 
duty  cycle  distortion.  Limiting  the  power  couple!  into  tne 
receiver  from  the  cable  will  prevent  this  from  securing. 

The  minimum  cable  length  may  be  determined  for  a 
specific  cable  by  calculating  the  power  launched  into  it  'ey 
usiag  the  typical  curve  from  Figure  5.4.  The  typical  value 
of  launched  power  should  be  used  since  this  will  give  the 
most  power  into  the  receiver,  resulting  in  the  desired 
minimum  cable  length.  Continuing  the  example  of  the  type 
5083  cable,  the  typical  launched  power  is  0.9  microwatts. 
The  minimum  loss  margin  is  found  using  the  typical  power 
launched  into  the  cable  and  the  maximum  input  to  :e 
receiver  which  is  1  microwatt  for  the  3713R  as  stated  in  tr.a 


electrical  specifications  for  the  receiver. 

LM  min  =  lOlog  (Pl/Pin,aax»  (11» 

LM  min  =  lOlog  (0.9  microwatt/1  microwatt)  =  -0.u6d3  (121 

The  minimum  cable  length  for  the  specified  caole  in  this 
example  is: 

Xmin  =  LM,min/L^  f  » 3 > 

Xmin  =  -0.46d8/7x10'3  dBm  =  -65.7m  (14) 


If  this  minimum  cable  length  is  too  long  for  a  particular 
application,  use  a  smaller  core  diameter  canle,  a  cable  with 
a  higher  loss,  or  the  power  launched  into  the  cable  may  be 
reduced  as  described  in  section  Il.i.I.c  by  adjusting  the 
transmitter  power  adjust.  The  negative  sign  in  the  minimum 
length  indicates  there  is  no  minimum  iength  limitation  or. 
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this  fiber  optic  cable  whan  used  in  con  junction. 
3713T/R  set. 


4  at  n. 


P.  TYPICAL  APPLICATIONS 

The  3713T/R  fiber  optic  link  solves  such  data  transmis- 
sion  problems  as  crosstalk,  ringing,  and  echoes. 
Electromagnetic  radiation  interference  is  avoided  when  using 
a  fiber  optic  data  link  in  high  noise  environments. 
Lightning  damage  to  cables  and  connected  equipment  can  be 
eliminated  where  fiber  optic  cables  replace  metallic  conduc¬ 
tors.  Fiber  optic  cables  are  not  sibject  to  shorting  out 
when  emersed  in  water. 
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appendix  S 
test  PROGRAM 

//JOES  1MB  JOB  {15  92, 0129)  , '  F  IS  HER  SMC.  1  399*  ,  CLASSIC 
//♦MAIN  ORG=NPGVM  1.  1  5  92P,  LIN  ES  =  ( 5 0) 

//♦FORMAT  PR,DDNAME=PLOT.  SYS 7ECr R , DES P  =  LOCAL 

//  EXEC  FKTXCLGP, FARM. LKED=' LIST, MAP, XREF' , REGION. GO=27COK 

//Si SOT1  DD  ON IT=  SYSDA, 5  PACE=  (CYL, (3,8) ) 

//SYS0T2  DD  ONII=SYSDA,5?  ACE=  (OIL,  (3, 8)  ) 

//Si  SLIN  DD  SPACE*  (6080,  (80,80)  )  ,UMIT  =  SYSDA 
//FORT.SYSIN  DD  D  SN=MSS.  S  Y5  3  .  NON  IMSL.  3  0  0  RCE  (FO  ORT)  ,DIS?=SKR 
//  DD  ♦ 

C  A  LARGE  NUMBER  OF  ARRAYS  ARE  UTILIZED 

C  WITH  THE  INTENTION  OF  MAINTAINING  EACH 

C  BIT  OF  ORIGINAL  DATA  OR  CALCULATION, 

C  READY  FOR  RECALL  AT  ANY  TIME  FOR  FURTHER 

C  COMPUT  AITON.  EQUIVALENT  ARRAYS  ARE  EMPLOYED 

C  3ECAUS  E  THE  PLOTTING  ROUTINE  *DRAWP'  IS 

C  NOT  ABLE  TO  HANDLE  COMPLEX  ARRAYS. 

C  ARRAYS  * IT3',,RTB,r ' ALA8*  AND  ’TITLE*  ARE 

C  USED  TO  GENERATE  THE  VERSATEC  PLOTTER  OUTPUT. 

C  OCEAN  ARRAYS  ARE  INDICATED  VIA  THE  SUFFIX 

C  'O',  AND  LAND  ARRAYS  ARE  INDICATED  WITH  'L'. 

COKPLEX*8  XXO(8192),YYO(3192)  , 

C  CSO  (8192)  ,  COO  (8192)  ,  S  XO  (3  192)  , 

C  CRO  (8  192)  ,  CTO  (8192)  ,  CUO  (3  1  92)  , 

C  C70  (8192)  ,CWO  (8192)  ,CPO(9192)  , 

C  Z  YO  (8  1 9  2)  ,  CO  MOL  (8  1  9  2)  ,COL(819  2i  , 

C  CUL  (8192)  ,  CVL  (8192)  ,  C  WL  (3  192)  , 

C  COPOL  (8192),XXL(8192t  ,YYL(8192>.  ,Z XL  (819  2)  , 

C  ZYL  (8192)  ,  CSL  (8192)  ,  3PO  (3  192)  , 

C  BMO  (8  19  2)  ,  BP L  (8192)  ,3ML(3192)  , 
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C  SPOL  (81  92 )  ,SMOL  (8  1 92)  ,SP3(8192|  , 

C  S PL  (8192)  ,  SMO  (8192)  ,CTL(3192l  , 

C  CPL  (8192)  ,  CRL  (8192)  ,  SML  (3  192) 

DIMENSION  TIME  (8192)  ,  FREQ  (8192i  , 

C  K0RK(16384) ,FRQ2(3192 ) 

INTEGERS  ITS  (12) /12*0/ 

REAL*4  RTB  (28)  /23*  0.  0/, 

C  CROO(  16364)  ,CPCO(16334)  ,11)03(1  5  334)  , 

C  0000(16384), CTOO(1  633  4)  ,0700  (15364)  , 

C  CWOO(  16384)  ,CTLL(1o334)  , 

C  CPLL(16  384),CRLL(163  3  4)  ,OOLL(15  384)  , 

0  CULL(  16384)  ,C7LL(1633  4)  , 

C  CWLL(16384)  ,  COMOLL  (15  384)  , C0P3L L ( 16384)  , 

C  SPOO  ( 1638 4  )  ,  SP LL(  1  6 3  3  4)  , 

C  S MOO (16  38  4),SMLL(16334) 

REAL  ALAS  (22)  /'PDXO»  ,  •  PDSf  0  1  ,  »  ST  1  0  *  , 

C  *ST20»,  *ST30»  ,  'STOO*,  »PHA3*, 

'C  'GAMO’^PDXL' ,  *PDYL»  ,  •  ST1  L  •  ,  •  S  T  2  L  »  , 

C  *  ST3L ' , *  S  T  OL' , *  PHAL’ ,  ' G AML  * , 

C  '  SPLO*  ,  *5?  LL'  ,  *  SMOO'  ,  »  S MLL •  , •  0 P 0 L  •  ,'CMOL'/ 

REAL*8  TIT  LE ( 1 2) 

EQOIV  ALENCS(TITLE  ( 1)  ,  RTB  (5)  )  ,  (0  930(1)  , 

C  CRO  ( 1 )  )  ,(CPOO(1),CPO(1)), 

c  (cooo  (i)  ,coo(i) ) ,  (croo(D  ,cro  (i>  i  ,  (cuoo  (i»  ,cuo  (i) ) 

C  (CWOO  (1 )  ,  CTO  (1 )),  (C7  30(1)  ,CV0(1)  )  ,  (CTLL(1)  ,CTL(1>) 

C  (CPLL  ( 1 )  ,C?L  (1)  )  ,  (C3LL(1)  ,CRL(1)  )  ,  (C0LL(1>  ,C0L(1)  ) 

C  (COLL  (1 )  ,  C  fJL  ( 1 )  )  ,  (C7L  L  ( 1)  ,  C7  L  ( 1 )  )  ,  (CWLL  (1)  ,CWL(1)  ) 

C  (COPOLL  (1)  ,C3  POL  ( 1 )  )  ,  (CCMOLL(I)  , 

C  COMOL  (1)  )  ,  (SPOO  (1)  ,  5P0  (It  )  , 

C  (SPLL  (1)  ,  SPL  (1)  )  ,  (SMO  0(1)  ,SM0(1)  )  ,  (SMLL(1»  ,SML(1)  ) 
DATA  XX  0, Y  YO, CS0,C03,  ZXO, Z YO, 03 D , CTO , COO, C70, CWO, 

C  XXL,YYL,ZX I, Z YL, CS  L, B PO, BMO , BPL , 

C  BML, SPOL,  5M0L,S?0, SPL  , SMO, 

C  S  ML,CTL  ,03  L,C0L,C7LfC0PGL,C7L, 


C  COPOL,COMOL/278523*p.  0,3.0)/ 

DATA  TIMS, FREQ,FRQ2/24575*0.0/ 

IFRAME=81 92 
NR=30 

FNR=FLOAT  (NS) 

P=0.0 

DO  70  L 1= 1  ,NR 

THE  DO  LOOP  ENDING  WITH  STATEMENT  70  ENABLES 
THE  PROGRAM  TO  PROCESS  A  LARGE  AMOUNT  OF  DATA 
BY  REPEATING  THE  PROCESS  IN  BLOCKS. THE  DATA 
POINTS  FROM  EACH  RON  THROUGH  THE  DO  LOOP  ARE 
ADDED  TOGETHER  AND  EVENTUALLY  AVERAGED  BY  THE 
NUMBER  OF  RUNS  THROUGH  THE  DO  LOOP.  'NR' 
REPRESENTS  THE  NUMBER  OF  DATA  SEQUENCES  TO 
BE  AVERAGED.  1  SEQUENCE  CURRENTLY  EQUALS 
8192  DATA  POINTS  FOR  EACH  CHANNEL  OR  256 
SECONDS  OF  DATA.  FOR  THIS  TEST  SINE  AND 
COSINE  HAVES  OF  *  HZ  HIL  BE  USED. 

DO  60  JJ=1,IFRAME 

XXO  (JJ)  =SI N(P*.785393  163| 

YYO  (JJ)  =CQS(P*. 785393  16 3i 
XXL  (JJ)  =SIN(P*.  785 393  16  3| 

YYL  (JJ)  =COS(P*.785  39  3  163* 

?  =  P  >  1 
CONTINUE 

THE  FOLLOWING  SECTION  GENERATES  THE  TIME 
AND  FREQUENCY  ARRAYS  AND  NORMALIZES 
THE  INPUT  PCM  DATA  TO  VOLTAGE  FORM  IN 
PREPARATION  FOR  THE  FAST  FOURIER  TRANSFORM 
TO  THE  FREQUENCY  DOMAIN. 

N=8 19  2 
FN=FLOAT (N ) 

DELTAT=1./32. 

T=FN*DELTAT 


DELTAF=1./T 
DO  20  J=1  ,  N 

TIME  (  J)  =DELTAT*FLOAT  (  J) 

FREQ(  J)  =DELTAF*FL3Ar  (  J) 

FEQ2  (J)  =ALOG10  (FREQ(J)  ) 

CONTINUE 

THE  NEXT  FOUR  STATEMENTS  PERFORM 
AN  FFT  ON  THE  INPUT  TIME  SERIES 
DATA. 

CALL  FOUST  {XXO,N,1  ,-1  ,0,WORK) 

CALL  FOUST  (YYO,N,  1  ,-1  ,0,WORK> 

CALL  FOUST  (XXL,N,1,-1  ,0,WORK) 

CALL  FOUST  (YY L,  N,  1  , -1  ,  0  ,  W  0  RKt 

DO  40  K4- 1  ,N 

XXO (KU)  =XX  0(K4) /FN 

YYO  (K 4)  =YYO(K4)  /FN 

XXL  (K4)  =XXL(K4)  /FN 

YYL  (K4)  =YYL(K4)  /FN 

CONTINUE 

THE  NEXT  LOOP  CALCULATES  THE 
’JNN  OR  MALI  ZED  SPECTRAL  DENSITIES 
FOR  SIN3LE  SITS  ORTHOGONAL 
COMPONENTS  OF  THE  GEOMAGNETIC  FIELD, 
THE  INDIVIDUAL  SITE  CROS  SPECTRA 
3 ETWEE  N  COMPONENTS,  THE  INDIVIDUAL 
ORTHOGONAL  COMPONENTS  OF  SIGHT  AND 
LEFT  CIRCULARLY  POLARIZED  SPECTRA, 
AND  TWO  SITE  CROSS  SPECTRA  FOR  RIGHT 
AND  LEFT  CIRCULAR  POLARISATION. 

DO  30  11=1  ,N 

ZXO  (II)  =ZXG(II)  ♦  (XXO  (  II)  *  CONJG(XXO  (II)  )  ) 
ZYO  (II)  =ZYO(II)  ♦  (YYO  (  II)  *CONJG  ( YYO  (II)  )  ) 
CSO  (II)  =CSC(II)  ♦  (XXO  (II)  *CONJ3  (YYO  (II)  )  ) 
ZXL(II)  =ZXL(II)  ♦  (XXL  (  II)  *  CON J3 ( XX L  (II)  )  ) 


ZYL  (II)  =ZYL(II)  ♦  (YYL  (  II)  *CQNJ3(YTL  (II)  )  ) 
CSL  (II)  =CSL(II)  ♦  (XXL  (II)  *C0NJ3  (YYL  (II)  )  ) 
3P0(II)  =XXO(II)  +  ((0.#  1.)*YY0(II)  ) 

BMO  (II)  =XXO(II)  -  (  (0.  ,  1.)  *YYO(II|  ) 

B?L  (II)  =XXL(II)  +  (  (0.  ,  1.)  *  YYL  (II)  > 

2ML  (I I)  =XXL(II)  -(  (0.  r  1.)* TIL  (II)  ) 

SPOL(II)  =S  FOL  (II)  ♦  (3PO(II)  *C0N  J  3 (BPL  (II)  )  > 
S MOL (II)  =  S  MOL  (II)  *■  (BMO  (II)  *C3N J3  (BML  (II)  )  » 
SPC  (II)  =SPO(II)  ♦  (BPO  ( II)  *  CON  j  3  (BPO  (II)  )  ) 
SPL(II)  =SPL(II)  +  (BPL(II)  *CONJ3(3PL  (II)  )  ) 
SMO  (II)  =SaO(II)  +  (BMO  ( II) *CONJ3 (BMO  (II)  )  ) 
SML(II)  =S  M  L  (II)  +  (BML  (II)  *  CONJ3  (BML  (II)  )  ) 
CONTINUE 
CONTINUE 

THE  NEXT  LOOP  NORMALIZES  THE  ABOVE 
SPECTRA  AND  CALCULATES  POWER  SPECTRAL 
DENSITIES. 

DO  33  13=1  ,N 
ZXO  (13)  =ZX  C  (13 )  *T/FNR 
ZYO  (13)  =Z Y 0(1  3 )  *T/FN3 
CSO  (I  3)  =CSC(I3)  *T/FN R 
ZXL  (13)  =ZXL(I3)  *T/FNR 
ZYL  (13)  =ZYL(I3)  *T/FN3 
CSL  (13)  =CS  L(I3)  *T/FM 3 
SFOL  (13)  =S  POL  (13)  *T/FNP. 

SPO  (13)  =SP0(I3)  *T/FM3 
SPL  (13)  =S  P  L  (I  3  )  *T/FNa 
SMO  (13)  =SM0(I3)  *T/FN 3 
S  ML  (I  3)  =SML(I3)  *T/FSR 
SMOL(I3)=S  MOL  (13)  *T/F  NR 
CONTINUE 

THE  NEXT  LOOP  CALCULATES  STOKES 
PARAMETERS  0  THROUGH  3  OF  THE 
INDIVIDUAL  SITE  ORTHOGONAL 
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C  COMPONENTS  AND  THE  COHERENCE  OF 

C  THE  PLANAR  CIRCULAR  POLARIZATION 

C  PARAMETERS  BETWEEN  SITES. 

DO  44  14=1  ,N 

CTC  (14)  =  (ZXO(I4)*ZYO(I4)|  *2./T 

CPC  (14)  =  (ZXO  (I4)-ZY0(I4))  *2./(T*CTO  (14)  ) 

CFO  (14)  =  (  4  *CSO  (14)  )  /  (  CTO  ( 14 )  *T) 

COC  (14)  =CSO(I4)  /(CS2RT  (ZXO  (14)  i  *CSQRT  (Z YD  (14)  )  ) 

CKO  (14)  =ATAN2  (AIMASCOO  (14))  ,  R5  AL  (COO  (14)  |  ) 

COO (14)  =CS  QRT  (COO (14 )  *CONJG  (COD  (14) ) ) 

CTO  (14)  =4.  34 29448* CLOG  (CTO  (14)  . 

CUC  (14)  =4.  342 9448* CLOG  (ZXO  (14)  t 

CVO  (14)  =4.  34 29448* CLOG  (ZXO  (14)  t 

CTL  (14)  =(ZXL(I4)*ZYL(I4)|  *2./T 

CPI  (14)  =(ZXL(I4)-ZYL(I4))  *2.  /  (T  *CTL  (14)  ) 

CEL  (14)  *  (4  *CS L  (14)  )  /(CTL(  14)  *T> 

CCL  (14)  =CSL(I4)  /(CS28T(ZXL(I4)  )  *CSQRT  (ZYL(I4)  )  ) 

CTL  (14)  =4.  3429448*CL0G  (CTL  (14)  ) 

CUL  (14)  =4.  3429448*CL0G  (ZXL  (I'4)  * 

CVL  (14)  =4.  3429448*CL0G(ZIL  (14)  » 

CD  L  (14)  =AT  AN2  ( AIMAG(COL  (14)  )  ,  REAL  (COL  (14)  )  ) 

COL  (14)  =CS  QRT  (COL  (14)  *C0NJ3  (COL  (14)  )  ) 

CQPOL  (14)  =  SPOL  (14)  /(C  SQRT  (SPO  (I  4)  )  *CSQBT  (5  PL  (14)  )  ) 
CCI-OL  (14)  =  CS  QRT  (COPOL  (14)  *CONJ3  (COPOL  (14)  )  ) 

CO  MOL  (14)  =SMOL  (14)  /(C  SQRT  (SHO  (14)  )  *CSQRT  (5  ML  (14)  )  ) 
CO  MOL  (14)  =  CS  QRT  (COMOL  (14)  *CONJ3  (COMOL  (14)  )  ) 

SFO  (14)  =4.  3429  448*CL0G  (SPO  (14)  l 
SPL  (14)  =4.  3429448*CL0G  (SPL  (14)  | 

SMO  (14)  =4  .  3429  448*CL0G(S3O  (14)  » 

SML  (14)  =4.  3429448*CL0G(S3L  (14)  i 
44  CONTINUE 

C  VERSATEC  PLOT  OP  CALCULATED  QUANTITIES 

NPTS=10./D  ELTAF+1 . 

C  *  NP7S '  DETERMINES  NUMBER  0?  POINTS 
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NECESSARY  TO  PLOT  THE  0  13  13  HEPTZ  RAN 
ITE  (2 )  =0 
ITE (3) =20 
ITE  (4)  =10 
ITB  (6)  =1 
ITE  (7)  =1 
ITB  (12)  =1 
ETB  (1 )  =  0 
RTE  (2)  =0 
RTB  (3 )  =  ALA  E  (1 ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAWP  (NPTS  ,  FRQ2  ,  COOD  ,ITB,  3TB) 

IT3  (6)  =1 
RTB  (3  )  =ALA  E  (2) 

READ  (5,  3000)  TITLE 

CALL  DRAW?  (NPT  S,  FRQ2 ,  CVOD , ITB, 3 TB) 

ITB  (6)  =1 
RTB  (3)  -  ALA  E) (3) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW  p  (NPTS,  FRED,  CPOD  ,  ITB,  3  TB) 

ITE  (6)  =1 
RTB  (3)  =  ALA  3  (4) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPT  S,  FREQ,  CROD  ,  ITB,  3  TB) 

ITB  (6)  =1 
RTE  (3)  =ALA  8(5) 

READ  (5,  30  0  0)  TI  TLE 

CALL  DRAW?  (N?TS,PRE2,  CROD  (2)  ,IT3,RTB) 

ITB  (6)  =  1 
RTB  (3)  =ALA  B(6) 

READ  (5,3000)  TITLE 

CALL  DRAW?  (NPTS,FRQ2  ,  CTOD  , ITB,  3  TB) 

ITB  (6  )  =  1 
RTB  (3)  =ALAB(7) 


READ  <5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS, FREQ,  CWOD  ,  ITB,  3  TB) 

ITS  (6)  =1 
STB  (3)=ALAB(9) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS, FREQ,  COOD  , ITB,  3TB) 

ITB  (6)  =  1 
RTB  (3  )  =  ALA  3  (5  ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,?RQ2,  CULL  , ITB,  3TB) 

ITB  (6)  =  1 

RTB  (3  )  =  ALA  B  ( 1  0  ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?TS,FHQ2,  CVLL,ITB,  3TB) 

ITS  (6)  =1 

F.TB  (3  )  =  AL  A  B  ( 1  1) 

READ  (5,  3C0C)  TITLE 

CALL  DRAWP  (NETS,  FRE2,  C?LL  ,  ITB,  3  T  B) 

ITB  (6)  =1 

RTB  (3  )=  ALA  3  (1  2) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NETS,  7RE2,  CRLL  ,  ITB,  3  TB) 

IT5  (6)  =1 

RTB  (3  )  =  ALA  3  (1  3  ) 

E  EAD  (  5,  3000)  TITLE 

CALL  DRAW?  (NPTS,  FRE2,  CRLL(2)  , ITB,  RTB) 

ITB  (6)  =  1 

RTB  (3  )  =  ALA  B  (1  4 ) 

READ (5,  30  0  0)  TITLE 

CALL  DRAWP  (NPTS,  FRQ2  ,  CTLL  ,  ITB,  3  TB) 

ITB  (6)  =1 

RTB  (3 )  3  ALA  B  (1  5  ) 

READ (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,  FREQ,  CWLL  ,  IT3,  3  1 3) 
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ITB  (6)  =  1 

ET3  (3)=ALA3(16) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,FRE2,  COLL  ,  1 TB,  3  TB) 

ITB  (6)  =  1 

RTB  <3)=ALAB(17) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,FRQ>,  SP  00  ,  ITB,  3  TB) 

ITB  (6)  =  1 

RTB  (3  )  =  ALA  3(18) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS  #  FRQ2  ,  SPLL,  3TB,  3TB) 

ITB  (6)  =1 

RTB  (3  )  =  ALA  B  (1  9  ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAWP  (NPTS,?RQ2,  3X00 , ITB, 3 TB) 

ITB  (6)  =1 

RTB  (3 )  =  AL  A  B  (20 ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,  FRQ2,  SELL,  ITB,  3TB) 

ITB  (6)  =  1 

RTE  (3  )  =  AL  A  B  (2  1 ) 

READ  (5,  30  0  0)  TI  TIE 

CALL  DRAW  ?  (NPTS,  ERE,},  COPOLL,  IT  3,  RTB) 

ITB  (6)  =  1 

RTB  (3)  =ALAB(22) 

READ  (5,30  0  0)  TITLE 

CALL  DRAWP  (NPTS,FRS2,  COMOLL,  ITB ,  RTB) 
3000  FORMAT (6A8 ) 

STOP 

END 

/* 

//GO.SYSIN  DD  * 

TEST  POWER  SPECTRAL  DENSITY  OF  SIN  '4  32 


IN  XXO, 30  AVGS,  (REF  DB  VS  LOG  FREQ) 

T2Sr  POWER  SPECTRAL  DENSITY  OF  EDS  4  112 
IN  Y  YO  ,30  AVGS  ,  (R  EF  DB  VS  LOG  FREQ* 

TEST  STOKES  1/STOKES  0  OF  SIN  432 
IN  XXO  AND  COS  4HZ  IN  YYO,3D  AVGS 
TEST  STOKES  2/STOKES  0  OF  SIN  43  2 
IN  XXO  AND  COS  4  H  Z  IN  YYO,33  AV35 
TEST  STOKES  3/STOKES  0  OF  SIN  43  2 
IN  XXO  AND  COS  4  H  Z  IN  YY0,3D  AV3S 
TEST  STOKES  0  OF  SIN  4 HZ  IN  XXO  AND 
COS  4HZ  IN  YYO  ,3  0  AVGS ,  (REF  DB  VS  LDil 
TEST  PHASE  OF  SIN  4H2  IN  XXD 
AND  COS  4HZ  IN  YYO,30  AVGS 
TEST  COHER  OF  SIN  4 HZ  IN  XXO 
AND  COS  4HZ  IN  YYO,30  AVGS 
TEST  POWER  SPECTRAL  DENSITY  OF  SIN  4  32 
IN  XXL, 30  A  VGS ,  ( R  EF  D3  VS  LOG  FREQ) 

TEST  POWER  SPECTRAL  DENSITY  OF  CDS  4  32 
IN  Y  YL  ,  30  A  VGS.  (  R  EF  D5  VS  LOG  FREQ) 

TEST  STOKES  1/SIOKES  0  OF  SIN  4HZ 
IN  XXL  AND  COS  4HZ  IN  YYL,30  A  VGS 
TEST  STOKES  2/STOKES  0  OF  SIN  4HZ 
IN  XXL  AND  COS  4HZ  IN  YYL,30  AVGS 
TEST  STOKES  3/STOKES  0  OF  SIN  4HZ 
IN  XXL  AND  COS  4HZ  IN  YYL, 3D  AVGS 
TEST  STOKES  0  OF  SIN  4 HZ  IN  XXL  AND 
COS  4HZ  IN  Y YL ,  3 0  AVGS,  (REF  DB  VS  LOG) 

TEST  PHASE  OF  SIN  4HZ  IN  XXL 
AND  COS  4 HZ  IN  YYL,  30  AVGS 
TEST  COHER  OF  SIN  4HZ  IN  XXL 
AND  COS  4 HZ  IN  YYL,  30  AVGS 

TEST  RT  CIRC  POLARIZATION  PSD  OF  SIN  4  HZ  IN 
XXO  AND  COS  4HZ  IN  YYO,  (DB  VS  LOG),  3D  AVGS 
TEST  RT  CIRC  POLARIZATION  PSD  OF  SIN  4  HZ  IN 


XXL  AND  COS  4HZ  IN  YYL,  (D3  VS  LOG)  ,33  AVGS 
TEST  LEFT  CIRC  POLARIZATION  PSD  OF  313  4  32 
IN  XXO  AND  COS  4HZ  IN  YYO  (DB  VS  LOG).  3  0  AvGS 
TEST  LEFT  CIRC  POLARIZATION  PSD  OF  SIN  4  HZ 
IN  XXL  AND  COS  4  H  Z  IN  YYL  (DB  VS  LOG),  30  AVGS 
TEST  COHER  RT  CIRC  POLARISATION  OP  SIN  4 HZ 
IN  XXO, XXL  AND  COS  4  HZ  IN  YYO, YYL,  30  AVGS 
TEST  COHER  LEFT  CIRC  POLARIZATION  OF  SIN  4hZ 
IN  XXO, XXL  AND  COS  4 H Z  IN  YYO, YYL,  30  AVGS 


APPENDIX  g 

MASS  SID  RAGE  PROGRAM 

//JDEB1M8  J03  (1592,  0129)  ,' FISHER  SMC.  1399'  ,  CLA5:=F 
//♦MAIN  ORG=NPGVM  1.  15  92P 
//  EXEC  FOP.TXCL3 

//FORT.  SYS  IN  DD  D  SN=MSS.  S  YS3  .  NON  IMSL.  5  OURC2  (FO  UR 1)  ,  DI  S?=SHR 
//  DD  ♦ 

C  THE  ARRAY  'IN'  WILL  BE  USED  13 

C  RECEIVE  THE  DATA  PASSED  FROM 

C  THE  SUBROUTINE  'READ*  AND  THEN 

C  TRANSFERRED  TO  THE  APPROPRIATE 

C  XXX  OR  YYY  ARRAY. 

INTEGER*2  IN(16) 

COMPLEX* 8  XXX(8192)  ,YYY(3192) 

DIMENSION  FREQ  (8192)  ,TIME(8192)  ,  W  3  R  K  ( 1 6  3  8  -'■ ) 

DATA  XXX,YYY/16384*  (3.  0,0.0)  / 

DATA  TIMS,FREQ/16334*3.0/ 

C  THE  FOLLOWING  SECTION  READS  THE 

C  FIRST  SIX  SECONDS  OF  COMPUTER  TAPE 

C  AND  DISCARDS  THIS  DATA. 

ITL=  19  2 

DO  55  JJ=1  ,  ITL 

CALL  RD(  20,  IN, 200, IP.  SC  ,IRR) 

55  CONTINUE 
IFF AME=8 192 
NR=20 

FNR  =  FLOAT  (N  R) 

DO  70  LI  =1  ,  NR 


C 

c 

c 


THE  NEXT  LOOP  READS  THE  COMPUTER 
TAPE  USING  THE  PROGRAM  PROVIDED 
BY  MR.  TIM  STANTON  OF  NAVAL  POST- 


c 


GRADUATE  SCHOOL, 


DO  60  JJ  =  1 , IPR A  ME 

CALL  RD  (20  ,  IS,  1000,  IREC,IJR) 

XXX  (JJ)  =IN  (2) 

YYY  (JJ)  =  IN  (3) 

5  0  CONTINUE 
N=S 192 
?N=FLOAT  (N) 

DELTAT=1  ./3  2. 

T=FN*D2LTA  T 
DELTAF=1 ./T 
DO  20  J=1,N 

TIME  ( J)  =  DELTAT  *FLOAT  ( J  ) 

FREQ  (J)  =DELTAF*FLOAT( J) 

C  THE  NEXT  4  STEPS  CONVERT  THE 

C  DATA  TO  VOLTAGE  AND  ENSURES  THAT 

C  NO  ERRONEOUS  DATA  HAS  BEEN  INTRODUCED 

C  INTO  THE  ARRAYS. 

XXX  (J)  =  (XXX  (J)  -2045.5)  *10./20'4  5.  5 
XXX  (J)  =REAL  (XXX  (J»  ) 

YYY  (J)  =  (YYY  (J)  -2045.5)  *10.  /2045.  5 
YYY  (J)  =RZAL  (YYY  (J)  ) 

20  CONTINUE 

CALL  FOURT  (  XXX, N, 1,-1,  0, WORK) 

CALL  FOUar  (YYY,  N,1,  -1,  0,  WORK) 

DO  40  K4  =  1,N 
XXX  (K4  )  =  XX  X  (K4 )  /FN 
YYY  (K4)  =  YY  Y  (K4)  /FN 
40  CONTINUE 

C  THE  NEXT  LOO?  APPLIES  THE  SYSTEM 

C  TRANSFER  FUNCTION  TO  THE  TRANSFORMED 

C  FREQUENCY  DOMAIN  DATA.  THE  TRANSFER 

C  FUNCTION  CONVERTS  VOLTS  TO  NANOTESLAS. 

DO  9  L=1,N 
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FRQ  =  FREQ  (L) 

1  IF  (FRQ.LE.  1  5.)  GO  ID  2 

XXX  (L)  =XXX  (L)  /( 105.5-3  .  14*FRQt 
YYY(L)  =YYY  (L)  /  ( 181. 32-  7.  53  3*FRQ) 

GO  TO  8 

2  IF  (FRQ.LE.  1  0.)  GO  TO  3 

XXX  (L)  =XXX  (L)  /  ( 5.958* PRQ-3  0.97 1 
YYY  (L)  =YYY  (L)  /  (7.  166*?  RQ-3  9 . 9 9» 

GO  TO  8 

3  IF  (FRQ.LE.  7. 5)  GO  TO  * 

XXX  (L)  =XXX  (L)  /  (3.49  2*?  RQ-5 . 31i 
YYY  (L)  =YYY  (  L)  / ( 4 . 25 2* F RQ-1 0 . 3 5 ) 

GO  TO  8 

4  IF  (FRQ.LE.  5  .)  GO  TO  5 

XXX  (L)  =  XXX  (L)  /  (2.6311  *  FRQ+Q.  14657) 

YYY  (L)  =YYY  (L)  /  (3 .0  1  2*P  RQ-1  .  55) 

GO  TO  3 

5  IF  (FRQ.LE.  3.)  GO  TO  5 

XXX  (L)  =XXX  (L) /  (2. 6311*  FFQ*  0.14557) 

7  YYY  (L)  =YYY  (I)  /  (2.732*FRQ) 

GO  TO  3 

6  XXX  (L)  =XXX  (L) /(2.72*FRQ) 

GC  TO  7 

8  CONTINUE 

9  CONTINUE 

THE  NEXT  WRITE  STATEMENTS  SEND 
THE  CONVERTED  DATA  TO  MSS 
FOR  FUTURE  MANX?  J  LATI  ON  AND  RECALL. 
THIS  NEXT  SET  OF  STATEMENTS  ALSO 
PROVIDES  THE  USSR  A  DIAGNOSTIC. 
WRITS(21)  XXX 

WRITE  (6,  *)  XXX(1)  ,XXX(3  192| 

WRITE  ( 2 1 ) YYY 

WRITE  (6,*)  YYY(1)  ,YYY(3  19  2* 


70  CONTINUE 
ENDFILE  21 
STOP 
END 

SUBROUTINE  RD  (IUN,  10,  IRS,  IREC,  IRQ) 

C 

C  THIS  PROCEDURE  FJRNISHED  3Y  MR.  TIM  STANTON , 

C  DEPARTMENT  OF  OCEANOGRAPHY. 

C 

C  READ  DATA  FROM  ION,  ALLI3N  ,  CHECK  &  3ZTUR 

C 

C  IUN=TA?E  NUMBER,  EG  20 

C  I  G=I  NTEGE  R*2  ARRAY,  16  LONG, 

C  (VALUES  0-4095,  SUBTRACT  2048t  *5/2028.  GIVES  VOLT 

C  IRS*  NUMBER  OF  RESINCS  ALLOWED  (ERRORS) 

C  I REC=  COUNTER  OF  RECORDS  (FRAMES  OF  DATA) 

C  BLOCK  512  BITS,  32  BITS  =  RECORD 

C  800  3  PI  TAPE  UMLABLED 

C  Ir.Q=  NUMBER  OF  ACIJAL  RESINCS  (ERRORS) 

C 

c 

INTEGER  *  2  10  (16)  ,1?  (16) 

DATA  IRR  /0/ 

IF  (IREC.EQ.O)  IS  =  0 
I  ER=0 

20  FORMAT  (16A2) 

IF  (IS.  NE.  0)  GO  TO  50 
READ  (I UN,  20,  END=9  00 1  I? 

I REC= IR  EC+  1 
40  IS=IS+1 

IF  (IS.  LT.  17)  GO  TO  5  0 
READ  (I UN  ,  20,  END=9  03 )  I? 

IS=1 

IP  EC  =  IREC+- 1 
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ICH=IMASK(I?(IS)  ,3,0)  <■  1 
WHITE  (6,5  5)  ICH,  IS,  I  UN, IR EC 
FORMAT  (•  RES  YNCING  ICH,  IS  ,  IUN  ,  IHEC  *  ,  41 3  > 


IF  (ICH.  NE.  1)  GO  TO  *3 
DO  100  1  =  1  ,  16 
IO(I)  =1  SHI  FT  (IF  (IS)  ,4) 

ICH=IMASK(IP(IS)  ,3, 3)  «■  1 
IF  (ICH.  SQ.  I)  GO  TO  33 
IER=IEH+  1 

WRITE  (6,70)  IUN,  I R3C  ,1,  ICH,IER 
FORMAT  ('  UNIT*  ,13,  'HE  CO  HD', IS, 'CHAN  5  DATA  CH  ',214 
$  'ERRORS  '  ,  17) 

IS=IS+1 

IF  (IS.  LI.  17)  GO  TO  1  00 
READ  (IUN,  20,  END=903l  IP 
IS=1 

I?.EC=IHEC+  1 
CONTINUE 


IF  (IER.EQ  .0)  GO  TO  1  50 
IRR=I RR  + 1 

IF  (IRS. LT. IRS)  GO  TO  123 
WRITE  (6,  1  10) 

FORMAT  ('  1  STOPPED  IN  SUB  RD  BECAUSE 
C  OF  IRR.  GT.  ',  16  ,  «  AT  L  110'  ) 

I3Q=IR3 

STOP 

CONTINUE 

WRITE  (6,130)  IREC ,  I R  R 

FORMAT  ('RESYNC  AT  ?3  AMZ'  ,1b  ,  '  WITH  TOTAL  ERRORS', 17 
I  ER=0 
IBQ=IR3 
GO  TO  50 


153  CONTINUE 
RETURN 

903  WRITE  (6,9  10)  IUN,I3EC 

910  FORMAT  ('1  END  OF  UNIT  ',13,'  AT  REC  *,I7) 

STOP 

END 

FUNCTION  ISHIFT  <IN,N?LQ  = 

C  RETURNS  SHIFTED  VALUE  OF  1*2  WORD  IN 

C  -VE  L  EFT,  +V  E  RIO  HT  SHIFT 

C 

INTEGER  *  2  IN 
IP=IN 

IF  (I  ?.  LT.  0)  IP=IP+65  536 
IF  (NPLC.LT.O)  GO  TO  30 
I  SHIFT-IP  /  (2**IA3S  (S?  LC)  t 
RETURN 

30  ISHIFT-IP*  (2**IA35  (N?LC)) 

IF  (ISHIFT. 3T.  655351  ISHI FT«M 0 0 (ISHIFT , 65 3 36 ) 

RETURN 

END 

FUNCTION  I  MASK  (IN,I3L,I8R) 

C  MASK  1*2  WORD  IN  OUTSIDE  BITS  I3L  S  133 

C 

INTEGER  *  2  IN, 10 
IO=IN 

IF  (IBR.EQ  .3)  GO  10  3  0 
IT= ISHI FT  (IN,  I3R) 

IO  =  IT 

50  I ?= I SHI FT  ( 10,  I  3L~1  5-1  3F) 

10  =  1? 

IMASK=ISHIFT  (10,15-IBL) 

RETURN 

END 

/* 


//GD. FT2 1F00 1  DD  ONIT=3330V,  ASV3 P=PUB4 A, DISP* (HE*  ,  CAILS)  , 

//  DSN=&SS.S1592.  GHDM  A, 

//  DCB=(3ECFM=VB5#5LKSIZE=4J96,LRECL=4092l , 

//  SPA  CE  =  (CYL,  (8,  8)  ) 

//GD.FT20F001  DD  UNIT=343 3 -4 , VOL =S ER  =  3 S DT 1  A , DISP= (OLD , ? ASS)  , 
//  LABEL*  (1,NL,,E  N)  , 

//  DCB  =  <3ECFM  =  FB,  LR ECL  =  32 , 3 LKSIZE=5 1 2 , D E N  =  2) 
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APPENDIX  H 
MAIN  PROGRAM 

//J0EK1M3  JOB  (15  92, 0129)  ,'FISHER  SMC.  1  399  •  ,  CLASS*  ? 

//♦MAIN  0RG=NPG7  M  1.  1  5  92P,  LIN  ES  =  ( 6  0) 

//♦FORMAT  PH , D  DM  A  ME=P  LOT.  SIS VECT  ?  ,  DE S T  =  LOCA L 

//  EXEC  FRTXCLGP,PARM.LKED=‘  LISI ,MAP, XREF* , REGION. GO=2700K 

//SISUT1  DD  U  N  IT  =  S  Y  S  D  A  ,  3  ?  AC  E=  (CYL,  (8,8)  ) 

//SYSUT2  DD  ONIT=SYSDA, SPACE*  (CYL,  (3,8)  ) 

//SYSLIN  DD  SPACE*  (6030,  (  80,30)  )  ,UNIT  =  SYSDA 
//FDRT.SYSIN  DD  D  SN=MSS.  S  Y5  3  .  NONI MSL.  5  DO  RCE  ( FOUR3M  ,  DIS?*SH3 


DETAILED  DESCRIPTIONS  OF  THE 
VARIOUS  STEPS  AND  LOOPS  ARE 
CONTAINED  IN  APPENDICES  B  AND  C. 
I  NT  EGER  *2  IN  (16) 

COMPLEX*8  XXO  (8192)  ,IYO(3  1921  , 

C  CSC(8192)  , COO  (3192)  ,ZXO  (3192)  , 

C  CRC  (8192)  ,  CTO  (3192)  ,  CUO  (3  1  92(  , 

C  CVO  (819  2)  ,  CWO  (3192)  ,  3  PO  (3  19  2)  , 

C  Z  YO  (8  1 9  2)  ,  COMOL  (8192)  ,COL(3192t  , 

C  CUL  (819  2)  ,CVL  (8192)  ,CWL  (3192)  , 

C  COPOL  (3192  ), XXX  (81921  ,  YY?  (8192*  , 

C  XXL  (819  2)  ,  YYL  (8192)  ,ZXL  (3  19  2)  , 

C  ZYL  (8192)  ,  CSL  (3192)  ,  B  PO  (3  1  9 2»  , 

C  BMO  (8192)  ,  BPL  (8192)  ,  3  ML  (3  19  2)  , 

C  SPCL  (8192)  ,SMOL  (9192)  ,SPD  (8132)  , 

C  S?L(8192)  ,  SMO(8192)  ,CTL(3192|  , 

C  C?L(8192)  ,  CSL  ( 8192)  ,  S  ML  (3  1 92) 
DIMENSION  TIME  (819  2)  ,  FEED  (3192i  , 

C  WORK  (16  38  4) ,?RQ2(3192) 

I NT EGER *4  ITB  (12)  /12*  0/ 


R2AL*4 


RTB  (28)  /23*  0.  0/ , CROO  ( 16  3S4)  f 


C  CPOO  (16384)  ,CUOO(16334)  , 

C  COOO (16384),CTOO(16334)  ,2 VOO (15  334)  , 

C  CWCO  ( 1633  4 )  ,  CTLL(  1 63  3  4)  , 

C  CPLL  {  16  38  4),CRLL(16334)  ,:OLL(15334)  , 

C  COLL  ( 1638  4  ) ,  C VLL(  1 63 3  4)  , 

C  CWLL(  16  384) , COMOLL  (15  384)  ,  COP  DLL  ( 16  3  84)  , 

C  SPOO(  16384) ,S?LL(16334)  , 

C  SMCO  t'16384),SMLL(1633  4) 

REAL  AL  A3  (22)  /*  PDXO'  ,  '  PDY  0*  , ' ST  1 0  *  , 

C  ,ST20»,  'ST30*  ,  'SIOO* ,  •  PHAQ* , 

C  1 GAMC  '  ,  1  P  D  XL '  ,  '  ?D  Y  L'  ,  '  SI1  L  •  ,  •  S  2  2L  •  , 

C  »5T3L*,  'STOL*  ,  *?HAL«  ,  »GAML*  , 

C  '  SPLC>»,  *S  PIL'  ,  'SMOO'  ,  'SMLL*  ,  •  CPDL*  ,  '  CMOL'  / 

REAL*  8  TI T L2  ( 1  2) 

Z  QOIV  AL2NCE(TITLE(  1)  ,  RTB(5) ) ,  (CROO  (1)  , 

C  CRO  (1  )  )  ,  (C  POO  (  1)  ,  CPO  (  1)  )  , 

C  (COOC  (1)  ,COO(1)  ),  (CTOO(I)  ,CT0  (1)  )  ,  (COOO(1»  ,COO(1)  ) 

C  (CWOO  (1)  ,  CWO(1)  )  ,  (CVDO(1)  ,C70  (1)  )  ,  (CTLL  (1)  ,CTL  (1)  ) 

C  (CELL  (1 )  ,  C  PL  (  1 )  )  ,  (CRL  L  ( 1)  , CRL ( 1 ) )  ,  (COLL  (1)  rCOL  (1)  ) 

C  (COLL  ( 1 )  ,CUL(1))  ,  (  CVL  L  ( 1)  ,CVL(1)  )  ,  (CMLL  (1 )  ,  CWL  ( 1)  ) 

C  (COPOLL  (1)  ,COPOL(1)  I  ,  (COMOLL  ( 1 1  , 

C  COKOL  (1)  )  ,  (SPOO  (1)  ,3?0  (1)  )  , 

C  (SELL  (1 )  ,  S  PL  ( 1)  )  ,  (SMOO  (1)  ,SM0(1|  |  ,  (SMLL  ( 1 )  ,  SML  ( 1 )  ) 
DATA  XXO,Y YO,CSO,COO,  ZXO, ZYO,CRO, 

c  cro,coo,cvc,cwo,xxx,  f yy, 

C  XXL,YYL,ZXL,ZYL,CSL,  BPO,BMO,8PL, 

C  B  ML, S  POL,  S  MOL,  SPO,S?L,SMO, 

C  SML,CTL,CR L,  COL  ,CV  L,  C  OPGL  ,  C  WL  ,  C  0  P  OL, 

C  COMOL/2949  12*  (0.0,  0.  3  )  / 

DATA  TIME,  FREQ, FR22/24575  *0.0/ 

ITL=1 92 

DO  55  JJ=1  ,ITL 

CALL  RD  (20  , IN, 200,  IREC,IS  R) 


55  CONTINUE 

IFF.AM  E=31  9  2 
NR=20 

FNP=FLO AT  (NR) 

T  HE  FIRST  ST  ATEMENT  S  OF  LOOP  70 
RECALLS  DATA  PREVIOUSLY  STORED  IN 
THE  IBM  3033  IASS  STORAGE  SYSTEM. 
DO  70  L  1=1  ,NR 
READ  (21)  XXX 
READ  (21)  YY  Y 
DO  43  11=1  ,IFR  AME 
XXL  (II)  =XXX(II) 

YYL  (II)  =YY  Y(II) 

43  CONTINUE 

DO  60  J  J=  1  ,IFRAME 

CALL  3D  (20  ,IN,  1  000  ,  I R  EC  ,1  RR) 

XXO  (J  J)  =1  N  (2) 

YYC  (JJ)  =IN  (3) 

50  CONTINUE 
N  =  S  19  2 
F N=FLCAT ( N ) 

DELTA  T=1./32. 

T=FN*  DELI  AT 
D5LTAF=1./T 
DO  20  J=1 , N 

TIME  ( J)  =DELTAT*FLOAT  (  J) 

FREQ  ( J)  =D  ELI  AF  *FL0  AT  ( J) 

FRQ2  (J)  =ALCS10  (FP.EQ(J)  ) 

XXO  (J)  =  (XXO(J)  -204  5. 5  )  *5.  /204  5.  5 
XXO  (J)  =REAL(XX0  (J)  ) 

YYO  (J)  =  (YYO(J)  -  204  5.  3  )  *5.  /204  5.  5 
YYO(J)  =  RE  A  L  (YY  0  (J)  ) 

20  CONTINUE 

CALL  FOORT  (XXO,  N,1  ,-1  ,0,  WORK) 


CALL  FOUST  (YYO,N,  1,-1  #0#WO3Kl 

DO  40  K4= 1  ,N 

XXO  (X4)  =X  X  C (K  4 )  /FU 

YYO  (K4|  =YYO(X4)  /F H 

CONTINUE 

DO  9  L=1,N 

?RQ=F  REQ  ( L ) 

IF  (FR  Q.  LE.  15.  )  GO  TO  2 

XXO  (L)=XXO  (L)  /  (-6  1.35  +17.06S*FR2  -  0.5389*r  RQ*?RQ) 
YYG  (L)  =  YY  0  (L)  /  (-6  3.43  5  +1  1 . 7 1  0  5  *F RQ-0 . 2327 9 *FRQ*FF.Q) 
GO  TO  3 

IF  (FRQ.  LE.  10.  )  GO  TO  3 

XXO  (L)  =XXO  (L) /(-13.496  +  4.9329  *FRQ) 

YYC(L)  =YYO  (L)  /  (-6  3.43  5  +1  1 . 71  05  *FRQ  -0 . 232  79  *FFQ*FP.O) 
GO  TO  8 

IF  (FP.Q.  LE  .  7.  5)  GO  TO  4 

XXO  (L)  -XX  0  (L)  /  (3.  947*  FBQ-9.368) 

YYC  (L)  =  YYO  (L)  /(4.  > 95 * F3Q-  1  1 . 83 7 ) 

GO  TO  8 

I?  (FHQ.  LE.  5.)  50  TO  5 

XXC  (L)  =  XX  3  (L)  /(2.  3 105  *FEQ -0.5557) 

YYC  (L)  -  YYO  (L)  /  (2.  3 00 *  FRQ-.  50) 

GO  TO  3 

XXO  (L)  =  XXO  (L)  /  (2.  7  23*  FRQ) 

YYC  (L)  =  YYO  (L)  /  (2.  7  17*  FRQ) 

GO  TO  8 
CONTINUE 
CONTINUE 
DO  30  11=1  ,N 

ZXC(II)  =ZXO(II)  ♦  (XXO  (  II)  *C0NJ3  (XXO  (II)  )  ) 

ZYO(II)  =ZYO(II)  +  (YYO  ( II)  *CONJ3  (YYO  (II)  )  ) 

CSO  (II)  =CSO(II)+  (XXO  (II)  *CONJ3  (YYO  (II)  )  ) 

ZXL(II)  =Z  X  L(I I )  +  (XXL(II)*CONJ3(XXL(II)  )  ) 

ZYL(II)  =ZY  L(II)  ♦  (YYL(  II)  *C3NJ3(  YYL  (II)  )  ) 


CSL  (II)  =CSL(II)  ♦  (XXL  { II)  0 CON  J 3  (XXL  (II)  )  ) 

EPO  (II)  =XXG(II)  + (  (C.  ,  1.)*YYO(IC|  | 

BMC  (II)  =XXO(II)  -((3.  ,  1 . >  <•  YYO  (XI >  ) 

BPL(II)  =XXL(II)  ♦(  (0.  ,  1.)  *YYL(II)  ) 

BML  (II)  =XXL(II)  -  (  (  3.  ,  1.)  *  XYL  (II)  ) 

SPCL(II)  =S  POL  (II)  ♦  (BPO  (II)  *CDN J 3  (3PL  (II)  )  ) 
SMOL(II)  =S  MOL  (II)  ♦  (3.1  0  ( II )  *C0SJ3  (BML  (II)  >  ) 

SPC  (II)  =SPO(II)  +  (3?j  (II)  *CONJ3  (BPO  (II)  )  ) 

SPL  (II)  *S  ?  L(I  I)  +  (BPL(II)*CONJG(BPL  (II))  ) 

SMC  (II)  =S  M  0  (II )  +  (3 MO  (  II)  *C0NJ3  (3MO  (II)  )  ) 

SML  (II)  sSML(II)  ♦  (3  ML  (II)  *CONJ3  (BML  (II)  )  ) 

CONTINUE 

CONTINUE 

DO  33  13=1, N 

2XC  (13)  =Z  X  0  (I  3 )  *T/PM3 

ZYC  (I  3)  =ZY  C(I3)  *T/FN3 

CSC  (T  3)  =CS  0 (13  )  *r/?N3 

ZXL  (13)  =ZXL(I3)  *T/P3B 

Z  YL  (1 3)  =  Z  Y  L  ( 1 3  )  *r/FN3 

csl  (13)  *csl(i 3)  *r/?»a 

S ? C L  ( 1 3 )  =S  POL  (13)  *?/P  HS 

s?o  (13)  =s  ?c (i 3 )  *r/F’;a 

SPL  (13)  =S  PL  (I  3)  *T/?n 
SMC  (I  3)  =S  MC(I  3)  *T/F»a 
SML  (13)  =SML(I3)  *r/PS3 
SMCL(I3)  =S  MOL  (13)  *T/?  NP 
CONTINUE 
DO  44  1 4=  1  ,N 

CTC  (14)  =  (ZXO  (14)  +  ZY0  (14)  i  *2. /I 

CPO  (14)  =  (ZXO  (14)  -ZYO  (14)  (  *2.  /  ( T  *CT0  (1 4)  ) 

CSC  (14)  =  ( 4  *CS0  (14)  )  /(CTO(  14)  *11 

COO  (14)  =CSO(I4) /{CS2ST  (ZXO  (14)  i  *CSQRT  (ZYO  ( 1 4)  )  ) 
CSC  (14)  =  ATAN2  (AIM*  3  (0  00  (14)  )  ,  RE  A  L  (COO  (14)  )  ) 

COO  (14)  =CSQRI  (COO  (14 )  *COS  JG(OOD  (14)  )  ) 


CTO  (14)  =4.  34 29443* CL 3  G  (CTO  (14)1 
CUO  (14)  =4  .  34  2  9  4  4  8* CL 3 G  (7X0  (14)  i 
CVO  (14)  =4.  34  2  9  4  4  8*CLDG  (770  (14)  | 

CTL  (14)  =  (Z  XL (14)  *Z  YL  (14)  i  *2./T 

CPL  (14)  =  (2  XL  (14)  -Z7L  ( 14}  |  *2.  /  (  P*CTL  (14)  ) 

CRL  (1 4)  =  (4*CSL  (14)  )  /(CTL(I4)  *T) 

C0L(I4)sCSL(I4)/(CS3RT(SXL(I4)  *CSQRT  (Z7L  (14)  )  ) 

CTL  (14)  =4.  342  9  4  48*CL3G  (CTL  (14)  ) 

C  U  L  ( 1 4 )  =4.  342S448*CL0G  (ZXL  (14)  i 
CVL  (14)  =4.  3429448*CL0G  (ZYL  (14) 

C»L  (14)  «ATAH2  (AIUi3(:  OL  (14)  )  ,R3AL  (COL  (14)  )  ) 
col (14)  =cscar  (col(I4i  *con  ig(C3l  (14) ) ) 

COPOL  (14)  =  SPO  L  (14)  /(I  SQPT  (5P0  ( E  4 1  )  *CSQRT(SPL  (14)  )  ) 
C0P0L  (14)  =CS2ST  (COPOL  (14)  *C0NJ3  (COPOL  (14)  )  ) 

COMQL  (14)  =  5M0L  (14)  /(ZSQS  P  (SMD  (14)  )  *CSQST(SSL  (14)  )  ) 
COKOL  (14)  =CSQRT  (CO  MOL  (I4t  *C0NJ3  (COMOL  (14)  )  ) 

SPO  (I  4)  =4.  342 9448* CLOG  (S?0  (14)  i 
SPL  (14)  =4  .  34  2  9  4  48* CLOG  (SPL  (14) 

SMC  (14)  =4  .  3429448*CL0G  (SMU  (14)  | 

S  ML  (14)  =4.  3429448*CL0G(31L  (14)  i 
CONTINUE 

N?TS= 1 0  .  /  D  ELI  A  F  +  1 . 

ITB  (2)  =0 
ITB  (3 )  =20 
ITB  (4)  =10 
ITB  (6)  =  1 
ITB  (7)  =1 
ITB  (12)  =1 
RTE  (1)  =0 
RTB  (2)  =0 
RTB  (3)  =  ALA  B(1) 

READ  (5,  3000)  TITLE 

CALL  DR  A  in' P  (NPTS,FRQ2f  CUOO  ,IT3rRT3) 

ITB  (6  )  =  1 


RTB  (3  )  =  ALA  3  (2 ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?TS,FR02,  CVCO  ,Ii’3f  .<  T') 

ITS  (6)  =1 
RT3  (3  )  =  AL  A  B  (3 ) 

READ  (5,  3000)  TITLE 

CALL  DRAW?  (NPTS,FRE2,  CPOO  ,ITB,  3TB) 

ITB  (6)  =1 
RTE  (3)  =ALA  3(4) 

READ  ( 5, 30  0  0)  TITLE 

CALL  DRAWP  (SPTS,  FRE2,  CROO  ,ITB,  RTE) 

ITB  (6  )  =  1 

RT3  ( 3)  =ALAB(5) 

READ  (5  #  30  0  0)  TITLE 

CALL  DRAWP  (RPTS,FRE2,  CROD  (2)  ,  I  T  3,  RTB) 

ITB  (6)  =  1 

RTB  (3)  =AL A  B(6) 

READ  (5,30  0  0)  TITLE 

CALL  DRAWP  (NPTS,FRQ2  ,  CTCD  ,IT3,»r3) 

ITB  (6)  =1 
RTE  (3)  =ALA  B(7) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAW?  (NPTS,FRE2,  CWCD  ,ITB,SI3) 

ITB  (5)  =  1 
RTB  (3)  =ALA  B(8) 

READ  (5 , 30  0  0)  TITLE 

CALL  DRAW?  (NPTS,FR22,  CCOD  ,I?3rRr3) 

ITB  (6)  =  1 
RTB  (3)  =AL  A  B{9) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAWP  (RPTS,FRQ2,  COLL, ITB, STB) 

ITB  (6)  =  1 

RTB (3) =ALA  3(10) 

READ  (5, 30  0  0)  TITLE 
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CALL  DRAWP  (N?TS,FRQ2,  C7LL  ,1TB,  3  TB) 

ITS  (6)  =  1 

RTB  (3)  =ALAB(1  1) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAWP  (NPTS,FRE2,  CPLL  ,113,  S  I  3) 

ITB  (6)  =  1 

RTB  (3)  =AL  A  3 ( 1  2 ) 

READ  (5,30  0  0)  TITLE 

CALL  DRAW?  (NPTS,FRE0,  CRLL  ,ITB, RT3) 

ITB  (6 )  =  1 

STB  (3)  =AL A  B ( 1  3 ) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAWP  (NPTS  ,  FR E2  ,  CE  LL  ( 2)  ,  I  T  3 ,  RTB) 

ITB  (6)  =1 

RTB  (3)  =ALA  B(1  4) 

READ  (  5f  30  0  0)  TITLE 

CALL  DRAWP  (NPTS,  FRQ2,  CTLL  ,  ITB,  3  2  S) 

ITS  (6)  *  1 

RTE  (3 )  =AL  A  B(1  5 ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (NPTS,  FREQ,  CWLL,ITS,  RTB) 

ITE  (6)  =1 

RTE  (3 )  =  AL A  E (1  6  ) 

READ  (5,  30 0  0)  TITLE 

CALL  DRAWP  (NPTS,  FRE2,  COLL  ,  ITB,  R  TB) 

IT3  (6)  =1 

RTB  (3)  =ALAB(17) 

READ  (5,  30 0  0)  TITLE 

CALL  DRAW?  (NPTS,  FRQ2,  SPOD , ITB, R I B) 

ITB  (6)  =1 

RTB  (3)=ALAB(18) 

READ  (5, 30  0  0)  TITLE 

CALL  DRAWP  (N?TS,FRQ2,  SELL,  ITB,  RTB) 

ITB  (6)  =1 


RTB  (3)  *AL  A  B ( 1  9 ) 

READ  ( 5,  3000)  TITLE 

CALL  DRAW P  (NPTS,FRQ2,  SMOO, ITS,  RTB) 

ITB  (6)  =1 

RTB  (3 )  =  AL  A  B(20 ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAWP  (NPTS,?RQ2,  SMLL,IT3,  3  i  3) 

ITB  (6)  =1 

RTB  (3)  =ALA3(2  1 ) 

READ  (5,  30  0  0)  TITLE 

CALL  DRAW?  (N?TS,FRE0,  COPOLL,  ITB,  P.TB) 

ITB  (6  )  =  1 

RTB  (3)  =ALA  3(22) 

READ  (5, 300  0)  TITLE 

CALL  DRAW?  (NPTS,FRE2,  CO  MOLL,  ITB  ,  RTB) 

3000  FORMAT  (6 A8  ) 

STOP 

END 

SUBROUTINE  RD  (IUN,  10,  IRS,  IREC,IRQ) 

C 

C  THIS  PROCEDURE  FURNISHED  81  HR.  TIN  ST  ANTClt , 

C  DEPARTMENT  OF  OCEANOGRAPHY. 

C 

C  READ  DATA  FROM  ION,  ALLIGN  ,  CHECK  £  RETURN 

C 

C  IUN=TAPE  NUMBER,  EG  20 

C  IG=INTE3E  R+2  ARRAY,  16  LONG, 

C  (VALUES  0-4095,  SUBTRACT  2  04  3 1  *5/2028.  GIVES  VOLIAGE 

C  IPS*  NUMBER  OF  RESINCS  ALLOWED  (ERRORS) 

C  IREC=  COUNTER  OF  RECORDS  (FRAMES  OF  DATA) 

C  BLOCK  512  BITS,  32  BUS  =  RECORD 

C  800  BPI  TAPE  UNLABLED 

C  IRQ*  NUMBER  OF  ACT J A L  RESINCS  (ERRORS) 
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$ 


INTEGER  *  2  10  ( 16>  ,IP  (16)  = 

DATA  IRR  /0/ 

I?  (IREC.  EQ.O)  15  =  0 
IER  =0 

FORMAT  (1  6  A2) 

IF  (IS.  HE.  0)  GO  ID  50 
READ  (IUN,  20,  END=9Q0)  IP 
IEEC=IREC+  1 
IS=IS ♦ 1 

IF  (IS.  LT.  17)  GO  TO  5  0 
READ  (IUN,  20,  END=900|  IP 
IS=1 

irsc=ipec+ i 

ICH=IMASK  (I ?(IS)  ,3,0)  f  1 
WRITE  (6,55)  ICH,  IS,  IUN, IS  SC 
FORMAT  ('  RESYNCING  ICH, IS, IUN, IREC  ',413 

IF  (ICH. NS.  1)  GO  TO  40 
DO  100  1  =  1  ,  16 
IO  (I)  =1  3H I  FT  ( I  F  (IS)  ,4) 

ICH=IMASK  (I  P(IS )  ,3,0)  f  1 
IF  ( ICH.  SQ .  I)  GO  TO  3  0 
IEP=IER*1 

WRITE  (6,70)  IUN, IREC  ,I,ICH,ISR 
FORMAT  (' UNIT1  ,13, 'RECORD' ,15,' CHAN  5  DATA 
'EFEORS  ',17) 

I3=IS+1 

IF  (IS.  LT.  17)  GO  TO  1  00 
READ  (IUN  ,  20 ,  END=9  00)  I? 

IS=1 

IE  EC= IR  EC ♦  1 
CONTINUE 


IF  (IZR.EQ  .0)  GO  TO  1  50 
IRR=IRR*1 

IF  (IRR.LT.IRS)  33  T3  123 
WRITS  (6,110) 

113  FORMAT  {’ 1  STOPPED  IS  SUB  RD  BECAUSE 
C  OF  IRR.GT.  ',16  AT  L  1 1  0*  ) 

IRQ=I RR 
STOP 

UD  CONTINUE 

WRITE  (6,130)  IREC  ,  1 3  R 

133  FORMAT  (’RESYNC  AT  FR  AME'  ,16,' WITH  TOTAL  ERRORS', 17) 
IER=0 
I3Q=I RR 
GO  TO  50 
153  CONTINUE 
RETURN 

903  WRITE  (6,9  10)  IUN,I3SC 

913  FORMAT  ('1  END  OF  3  NIT  ’,13,’  AT  RSC  ’,17) 

STOP 

END 

FUNCTION  I  SHIFT  (IN,N?LC) 

C  RETURNS  SHIFTED  VALUE  OF  1*2  WORD  IN 

C  -VE  LEFT,  ♦V2  RI3HT  SHIFT 

C 

INTEGER  *  2  IN 
IP*  IN 

IF  (IP.LT.0)  IP=IP  +  55  536 
IF  (NPLC.LT.3)  GO  TO  30 
I3HIFT=IP/  (2**IAa  (N?  LC)  ) 

RETURN 

30  ISHIFT=IP*  (2**IABS  (NPLC) ) 

IF  (ISHIFT.GT. 65535)  ISHI FT=M0  D (ISHI FT, 65336) 

RETURN 

END 
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FUNCTION  I  MAS  K  (IN,IBL,ISR) 

C  MASK  1*2  WORD  IN  OUTSIDE  BITS  IBL  S  IE? 

C 

INTEGER  *  2  IN, 10 
IO  =  IN 

IF  (IBR.EQ.O)  GO  TO  50 
IT  =  ISHIFT  (IN,IBR) 

I  C= IT 

50  IP  =  I  SHI  FT  ( ID  ,  IBL-1  5- 1  BR) 

IO=IP 

IMASK=ISHIFT(I0,15-IBL) 

RETURN 

END 

/* 

//GO  .FT20F001  DD  UNIT=340 0- 4 , VOL =SER=G MDT 1 3 , DIS ?= ( OLD, P A 
//  LABEL*  (1,NL,,  IN)  , 

//  DC 3=(RECFM=F3,  LF.EC L  =  3 2 , 3  LKSIZE=5  1 2  ,  DEN*  2) 

//G0.FT21F001  DD  ONIT=3330V,  MSV3P=PU34 A, DISP* (OLD,  KEEP)  , 
//  DSN=MSS.S1592.  GMDT1  A, 

//  DCB  =  (RECFM  =  VB5 , BLK SIZE=4 0 96 , LRECL*4 0 9 2) 

//GO  .  S  YSDUMP  DD  SYSOUT=A 
//GO.SYSIN  DD  * 

X  COIL  PSD,  MTRY  3A7 , 1 122-1 2 54  LOCAL , 1 7AUG82 , 

20  A  VGS,  32S/SEC, 5VQLT,  (DB  NT**2  V5  LOG  FREQ) 

Y  COIL  PSD,  MTRY  BA Y,  1  122-1  2  54  LOCAL, 1 7A UG82 , 

20  A  VGS,  22S/SEC,  5V0LT,  (DB  NT**2  V3  LOG  FREQ) 

STOKES  NO.  ONE,  MTRY  BAY,  1  1  22-1  254 
LOCAL,  17AUG32,  2  0  AVGS,  32S/SEC,  5 VOLT 
STOKES  NO.  TWO,  MTRY  BAY,  11  22-1  254 
LOCAL,  17AUG82,  20  AVGS,  32S/SEC,  5VOLT 
STOKES  NO.  THREE,  MTRY  BAY,  1122-1254 
LOCAL,  1  7  AUG  82  ,  2  0  AVGS,  32  S  /SEC  ,  5VOLT 
STOKES  NO.  ZERO,  MTRY  BAY,  1  122-  1254  LOCAL, 

17A3G82,  20  AVGS,  32S/SEC ,  5  VOLT ,  (LOO  VS  LOG) 
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PHASE  OF  X  AND  Y ,  MTRY  BAY  1 122-1254 
LOCAL,  17AUG82,  2  0  AVGS,  325/SEC,  57DLT 
COHER  OF  X  AND  Y  COILS,  MTRY  BAY,  1122-1254 
LOCAL,  17AUG82,  2  0  AVGS,  32 S /SEC,  5VDLT 
X  COIL  PSD,  LA  MESA, 11 22-12545 LO CAL, 17A0G82, 

20  AVGS,  32  S/S  EC,  107  OLT,  (DB  NT**2  YS  LOG  FREQ) 

Y  COIL  PSD,  LA  MESA,  1  122-1254  LOCAL,  1  7  AUG82  , 

20  AVGS,  32S/SEC,  10VOLT,  (DB  NT*  *  2  7S  LOG  FREQ) 
STOKES  NO.  ONE,  LA  MESA,  1  122012  54 
LOCAL,  17AOG82,  2  0  AVGS,  325  /SEC,  10VDLT 
STOKES  NO.  TWO,  LA  MESA,  11  2  2-12  54 
LOCAL,  17AUG82,  20  AVGS,  325/SEC,  1370LT 
STOKES  NO.  THREE,  LA  MESA,  1  122-1254 
LOCAL,  17AOG82,  20  AVGS,  32 S /SEC ,  13Y0LT 
STOKES  NO.  ZERO,  LA  MESA,  1122-1254  LOCAL, 
17A3G82,  20  AVGS,  32S/SEC ,  10VOLT,  (LOG  VS  LOG) 
PHASE  OF  X  AND  Y,  LA  MESA,  1  122-1254 
LOCAL,  17AUG82,  20  AVGS,  325/SEC,  10VDLT 
COHER  OF  X  AND  Y  COILS,  LA  MESA,  1122-1254 
LOCAL,  17AUG82,  2  0  AVGS,  323  /SEC,  10YDLT 
RT  CIRC  POLARIZATION  PSD, MTRY  BA  Y  ,  1 1  22- 1  254 
LOCAL, 17AOG32, 20  AVGS ,32S/3 E C, 57 OLT, (D B  VS  LOG) 
RT  CIRC  POLARIZATION  PSD,  LA  MESA,  1122-1254 
LOCAL,  17ACG82, 20  AVG S , 32S /5 E C , 10 VOLT,  ( D B  VS  LOG) 
LEFT  CIRC  POLARIZATION  PSD,  MTRY  5 A Y ,  1  1  2 2- 1  2 54 
LOCAL, 17ADG82,  20  AVGS ,32  S/5  EC ,5  VOLT ,( DB  VS  LOG) 
LEFT  CIRC  POLARIZATION  PSD,  LA  MESA,  1  122-1254 
LOCAL, 17A0G82, 20  AVG S , 32S /5 S C, 10 VOLT, ( DB  VS  LOG) 
COHER  RT  CIRC  POLARIZATION  MTRY  BAY/LA  MESA, 
1122-1254  LOCAL,  17AUG82,  20  AVGS,  325/SEC 
COHER  LEFT  CIRC  POLARIZATION  MTRY  BAY/LA  MESA, 
1122-1254  LOCAL,  17AOG82,  20  AVGS,  325/SEC 
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